NASA/TM-2012-216001 



Shielded-twisted-pair cable model 
for chafe fault detection via 
time-domain reflectometry 


Stefan R. Schuet 
Ames Research Center 
Moffett Field, California 

Dogan A. Timugin 
Ames Research Center 
Moffett Field, California 

Kevin R. Wheeler 
Ames Research Center 
Moffett Field, California 


March 2012 


NASA STI Program ... in Profile 


Since its founding, NASA has been 
dedicated to the advancement of 
aeronautics and space science. The 
NASA scientific and technical 
information (STI) program plays a key 
part in helping NASA maintain this 
important role. 

The NASA STI Program operates 
under the auspices of the Agency Chief 
Information Officer. It collects, 
organizes, provides for archiving, and 
disseminates NASA’s STI. The NASA 
STI Program provides access to the 
NASA Aeronautics and Space 
Database and its public interface, the 
NASA Technical Report Server, thus 
providing one of the largest collection 
of aeronautical and space science STI 
in the world. Results are published in 
both non-NASA channels and by 
NASA in the NASA STI Report Series, 
which includes the following report 
types: 

• TECHNICAL PUBLICATION. 
Reports of completed research or a 
major significant phase of research 
that present the results of NASA 
programs and include extensive data 
or theoretical analysis. Includes 
compilations of significant scientific 
and technical data and information 
deemed to be of continuing reference 
value. NASA counterpart of 
peer-reviewed formal professional 
papers, but having less stringent 
limitations on manuscript length and 
extent of graphic presentations. 

• TECHNICAL MEMORANDUM. 
Scientific and technical findings that 
are preliminary or of specialized 
interest, e.g., quick release reports, 
working papers, and bibliographies 
that contain minimal annotation. 
Does not contain extensive analysis. 

• CONTRACTOR REPORT. 

Scientific and technical findings by 
NASA-sponsored contractors and 
grantees. 


• CONFERENCE PUBLICATION. 
Collected papers from scientific and 
technical conferences, symposia, 
seminars, or other meetings 
sponsored or co-sponsored by NASA. 

• SPECIAL PUBLICATION. 

Scientific, technical, or historical 
information from NASA programs, 
projects, and missions, often 
concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. 
English- language translations of 
foreign scientific and technical 
material pertinent to NASA’s 
mission. 

Specialized services also include 
creating custom thesauri, building 
customized databases, and organizing 
and publishing research results. 

For more information about the NASA 
STI Program, see the following: 

• Access the NASA STI program home 
page at http://www.sti.nasa.gov 

• E-mail your question via the Internet 
to help@sti.nasa.gov 

• Fax your question to the NASA STI 
Help Desk at 443-757-5803 

• Phone the NASA STI Help Desk at 
443-757-5802 

• Write to: 

NASA STI Help Desk 
NASA Center for AeroSpace 
Information 
7115 Standard Drive 
Hanover, MD 21076-1320 



NASA/TM-2012-216001 



Shielded-twisted-pair cable model 
for chafe fault detection via 
time-domain reflectometry 


Stefan R. Schuet 
Ames Research Center 
Moffett Field, California 

Dogan A. Timugin 
Ames Research Center 
Moffett Field, California 

Kevin R. Wheeler 
Ames Research Center 
Moffett Field, California 


National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett Field, California 94035-1000 


March 2012 


Acknowledgments 

This work was funded by NASA’s Aeronautics Research Mission Directorate under the 
Vehicle Systems Safety Technologies Project of the Aviation Safety Program. 


The use of trademarks or names of manufacturers in this report is for accurate reporting and 
does not constitute an offical endorsement, either expressed or implied, of such products or 
manufacturers by the National Aeronautics and Space Administration. 


Available from: 


NASA Center for AeroSpace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 
443-757-5802 






Contents 


List of Tables 3 

List of Figures 4 

List of Symbols 5 

Executive Summary 8 

1 Introduction 9 

1.1 Problem statement 9 

1.2 Literature survey 10 

1.3 Report overview 11 

2 Theory 13 

2.1 Basic model formulation 13 

2.1.1 Cable specification 13 

2.1.2 Maxwell equations 14 

2.1.3 Field solutions 16 

2.1.4 Eigenvalue equation 16 

2.2 Uniform cable 17 

2.2.1 Eigenvalues 17 

2.2.2 Mode fields 17 

2.2.3 Choice of r c 20 

2.2.4 Exact potentials 22 

2.3 Equivalent circuit model 23 

2.3.1 Basic variables 23 

2.3.2 Transmission-line equations 24 

2.3.3 Calculation of C 26 

2.3.4 Calculation of Z 28 

2.3.5 General solution 32 

2.4 Twisted cable 36 

2.4.1 Eigenvalues 36 

2.4.2 Mode fields 36 

3 Simulation 38 

3.1 Overview 38 

3.2 Geometric fidelity 38 

3.3 Mesh 39 

3.4 2D analysis 40 

3.5 3D analysis 40 

3.6 Chafe faults 43 


1 



4 Experiment 46 

4.1 Overview 46 

4.2 System model 46 

4.2.1 Common-mode cable model 46 

4.2.2 A simple chafe fault model 47 

4.2.3 System S'-paranreters 48 

4.3 Bayesian parameter inference 49 

4.4 Experimental results 51 

4.4.1 Nominal cable examples 51 

4.4.2 Faulty cable example 57 

5 Conclusion 60 

5.1 Summary of main results 60 

5.2 Areas of future research 61 

5.3 Acknowledgments 61 

5.4 Contact information 61 

Bibliography 62 


2 



List of Tables 


2.1 Shielded-pair cable parameters 14 

2.2 Comparison of SUP cable capacitance values obtained with different models 28 

3.1 Capacitances from 2D analysis 40 

3.2 Eigenvalues from 3D analysis 41 

3.3 STP cable eigenvalues from 3D analysis for theoretical model 41 

3.4 Capacitances with faults for the even mode 44 

4.1 Parameter prior information 51 

4.2 Parameter estimates for example 1 52 

4.3 Parameter estimates for example 2 54 


3 



List of Figures 


2.1 Shielded-pair waveguide cross-section 13 

2.2 A pair of helically wound current-carrying filaments inside a cylindrical shield 15 

2.3 Cut-off frequencies of the even and the odd waveguide modes of SUP cable 18 

2.4 Potential distributions for the even and the odd modes of SUP cable 19 

2.5 Surface charge distributions for the even and the odd modes of SUP cable 19 

2.6 Electric-field distributions for the even and the odd modes of SUP cable 21 

2.7 Optimal filament positions for the even and the odd modes of SUP cable 22 

2.8 Variation of the approximate potential on the periphery of conductor 1 for the even and the 

odd modes of SUP cable 22 

2.9 Comparison of the approximate analytical and the exact numerical solutions for the even and 

the odd mode potentials of SUP cable 23 

2.10 Magnitudes and phases of the surface impedances of a round wire and a cylindrical tube in a 

coaxial arrangement 30 

2.11 Real and imaginary parts of the even- and the odd-mode eigenvalues for lossy SUP cable ... 33 

2.12 Phase velocities and attenuation lengths of the even and the odd modes for lossy SUP cable . 33 

2.13 Magnitudes and phases of the even- and the odd-mode characteristic impedances for lossy 

SUP cable 34 

2.14 A transmission-line circuit featuring STP cable 34 

3.1 Cross-sections of simulated cable geometries 38 

3.2 Meshes for STP cable and hole fault 39 

3.3 Simulated electric-field distributions for the even mode 40 

3.4 Simulated eigenmode fields 42 

3.5 Electric field penetration into metals 42 

3.6 Attenuation per meter in theoretical STP cable 43 

3.7 Profile of transverse E x field 44 

3.8 Profile of transverse E y field 45 

3.9 Electric-field distribution for faulty cable 45 

3.10 Axial electric field in the hole 45 

4.1 Model of chafe fault on a shielded cable 47 

4.2 VNA hardware model 49 

4.3 .S'-parameter block diagram of the measurement setup 49 

4.4 Example 1: Samples from the marginal posterior distribution 53 

4.5 Example 1: Comparison between the measured and modeled real part of 5n 53 

4.6 Example 2: Samples from the marginal posterior distribution 55 

4.7 Example 2: Comparison between the measured and modeled real part of Sn 55 

4.8 Example 2: Comparison between the measured and modeled real part of S 21 56 

4.9 Comparison of modeled and measured attenuation constants 56 

4.10 Frequency-domain comparison between nominal and chafed models 57 

4.11 Input interrogation signal used for the time-domain analysis 58 

4.12 Time-domain comparison between nominal and chafed models 58 


4 



List of Symbols 


A 

A 

B 

c 

C 

e 

c 

d 

e 

E[-] 

E 

/ 

F 

9 

G 

G 

H 

h£\-) 

i 

i 

I 

3 

imO) 

j 

j 

J 

Jm(-) 

k 

K 

Km(-) 

l 

i 

L 

in 

M 

n 

n 

N 

A f + 

N m (-) 

P 

P(') 


unknown coefficient 
insulator cross-section 
unknown coefficient 

speed of light in the insulator; per-unit-length capacitance 

unknown coefficient 

closed contour 

capacitance matrix 

distance 

2.71828 

statistical expectation operator 
electric field 
(linear) frequency 

magnetic flux; normalized potential 
measurement forward model 
measurement gain 
conductance matrix 
magnetic held 

Hankel function of the first kind 

conductor current 
current vector 
constant current 
imaginary part 

modified Bessel function of the first kind 

conductor index 

current vector 

volume current density 

ordinary Bessel function of the first kind 

wavenumber in the insulator 

surface current density 

modified Bessel function of the second kind 

per-unit-length inductance 

length 

inductance matrix 

azimuthal eigenvalue; measurement index 

number of measurements 

Monte-Carlo sample index 

unit normal vector 

number of samples 

nonnegative Gaussian distribution 

ordinary Bessel function of the second kind 

helix period 

probability density function 


5 



y 

q 

Q 

r 

R 

5ft 

s 

Rmn 

s 

t 

u 

u 

v 

V 

V 
w 
W 

X 

y 

y 

z 

z 

Z 

( x,y,z ) 


integration path 

wire center offset from origin 

electric charge 

radius 

(effective) radius 
real part 

noise standard deviation 
S'-parameter 
S matrix 
time 

voltage vector 
mode amplitude 

conductor voltage; propagation velocity 
voltage vector 
mode amplitude 
width 

bounding wall 
eigenvector 
measurement vector 
eigenvector 

per-unit-length impedance 
unit axial vector 
impedance matrix 
rectangular coordinates 


a 

P 

1, 7 m 

r 

5 {-) 

$mn 

A(j> 

At 

Az 

e 

£ 

c 

r? 

6 

K 

A 

y 

V 

t 

7T 

n 
p , e 

a 

S 

T 

<t>, <P 
$ 

X 

tjj 


attenuation index 

phase index 

radial eigenvalue 

reflection coefficient 

Dirac delta 

Kronecker delta 

azimuthal interval 

time interval 

axial interval 

permittivity 

relative permittivity 

surface impedance 

impedance coefficient 

unknown parameter vector 

helix twist rate 

axial eigenvalue 

permeability 

measurement noise 

wavenumber in metal 

3.14159 

Hertz vector 

unit radial vector 

conductivity 

surface charge density 

capacitance coefficient 

unit azimuthal vector 

electrostatic potential 

eigenvalue 

helix pitch angle 


6 



& 

Debye potential 

UJ 

(angular) frequency 

(p, <t>,z) 

(global) cylindrical coordinates 

(p, V, z) 

(local) cylindrical coordinates 

V 

nabla 

II 

axial component 

_L 

transverse component 

AC 

alternating current 

DC 

direct current 

EWIS 

electrical wiring and interconnect system 

fft 

fast Fourier transform 

if ft 

inverse fast Fourier transform 

MCMC 

Markov-chain Monte Carlo 

RF 

radio frequency 

STP 

shielded twisted pair 

SUP 

shielded uniform pair 

TDR 

time-domain reflectometry 

TDT 

time-domain transmissometry 

TEM 

transverse electromagnetic 

VNA 

vector network analyzer 


7 



Executive Summary 


Problem statement. Studies show that electrical wires on aircraft often experience chafing, which causes 
short circuits or broken wires leading to loss of system functionality, accompanied by smoke and fire events, 
presenting a serious threat to flight safety. The capability to locate and characterize chafing can potentially 
enable preventive maintenance well before failures occur, thus maximizing safety while minimizing out-of- 
service time. In order to realize this capability, however, highly accurate and computationally efficient models 
are needed to detect, in noisy operating environments, the small electrical effects that chafing can cause. 

Overview. This report details the development, verification, and validation of an innovative physics- 
based model of electrical signal propagation through shielded-twisted-pair cable, which is commonly found 
on aircraft and offers an ideal proving ground for detection of small holes in a shield well before catastrophic 
damage occurs. The accuracy of this model is verified through numerical electromagnetic simulations using a 
commercially available software tool. The model is shown to be representative of more realistic (analytically 
intractable) cable configurations as well. A probabilistic framework is developed for validating the model 
accuracy with reflectometry data obtained from real aircraft-grade cables chafed in the laboratory. 

Main results. The analytical modeling process and its subsequent verification through rigorous numerical 
simulations and its validation through laboratory measurements resulted in a number of insights relevant 
to the practice of chafing fault detection. Firstly, the various assumptions and approximations made in 
deriving the theoretical model were found to be acceptable for the typical cable geometries. In particular, 
the numerical simulations have shown that the eigenvalues derived from the theoretical model are close 
enough to the actual eigenvalues of the impedance-controlled cable to be realistic for practical application. 
This demonstrated ability to derive a highly accurate yet closed-form physics model is a key contribution of 
this work, as it is crucial for the subsequent development of fast and reliable fault detection algorithms. 

The laboratory validation of model accuracy on real aircraft-grade cables required additional steps to 
model the measurement setup itself and to employ probabilistic methods to infer the key underlying param- 
eters such as effective dielectric permittivity of the insulation and the finite conductivity of the wires and 
the shield. Notably, the optimal estimates inferred from data differed from the manufacturer’s specifications. 
Overall, the analytical model was able to match lab measurements to within 2 % absolute maximum error 
out to 1 GHz. Not surprisingly, it was observed that the slight analytical correction for twisting does not 
have a significant effect on the ability to fit the model to measurements, though it did seem to reduce the 
bias on the parameter estimates. It was therefore judged that other (unmodeled) elements in the system 
( e.g ., cable imperfections, connectors, etc.) are mainly responsible for the residual. 

The analytical derivation and supporting simulations lay the foundation for a trade-space study to aid 
in the requirements specification of both a maintenance tool and an on-board fault diagnosis system. One 
important observation, relevant to any such study, is that for a fixed fault size and distance from the source, 
the angular fault location relative to the two inner conductors has a strong effect on the scattered field, and 
hence creates ambiguity in the ability to quantify fault location and size. It was confirmed that the simple 
impedance-discontinuity model of a chafe is able to replicate the tell-tale “differentiation” effect of a small 
hole on a propagating pulse. It was also established empirically that average-size chafes do not produce 
a sensible reflection or transmission signal below 200 MHz. All other classical trade-offs also exist. For 
example, a fast rise-time interrogation pulse will be better able to discern the edges of a small fault by using 
higher frequencies, with the down-side that high frequencies are rapidly attenuated in long cables. 



Chapter 1 


Introduction 


1.1 Problem statement 

The electrical wiring and interconnect system (ewis) in any vehicle is a critical, and sometimes overlooked, 
subsystem where relatively minor issues can grow and eventually lead to serious safety problems like smoke, 
fire, and loss of critical system functionality. As such, these types of faults are sometimes considered either 
single-point events where a single component {e.g., a wire) fault causes a system to fail, or common-cause 
events where a single fault causes multiple systems to fail. Such events are generally considered high risk 
no matter how improbable they may be. Wiring faults, however, do frequently occur in practice, and are 
clearly capable of leading to both single-point and common-cause fault events; for examples in aviation, see 
[ 1 , 2 ].* 

While the technology for detecting hard faults {e.g., opens, shorts, and arcing) is available, it only enables 
mitigation after serious sustained electrical issues occur. Preventive measures, on the other hand, require 
the detection of precursor problems. Wire chafing is one such precursor to both intermittent and full-system 
breakdown, and is among the most commonly occurring issues identified in wiring health maintenance 
studies [1]. Thus, the ability to detect and locate chafe faults enables the repair of wiring systems well before 
serious problems occur, and even in some cases when intermittent faults observed in operation disappear 
during maintenance. To date, the failed efforts targeted at detecting chafing faults in practice have focused 
primarily on hardware development. As a result, there is almost no existing research detailing the underlying 
physics needed to thoroughly understand how signals propagating through the interconnect system are 
affected by these faults. Such an understanding would allow one to design fault detection algorithms that can 
optimally identify and locate chafes by extracting fault signatures from native or injected signals propagating 
on the cable. 

This work seeks to fill the primary technology gap identified above by detailing a physics-based model for 
signal propagation in shielded-twisted-pair cable, commonly used in aircraft, among a great variety of other 
application domains. The chafing of the insulation on an unshielded cable cannot be detected remotely, 
since the permittivity contrast between most dielectrics and the ambient medium {e.g., air) is too small. 
Furthermore, the continuously changing distance between an unshielded cable and some nearby metallic 
surface that serves as a ground plane produces spurious signal reflections which are typically much greater 
than those caused by chafing alone, as was demonstrated in [3]. The shield not only isolates the internal 
signals flowing down the cable from the noisy environment, but also serves as a uniform conductor for the 
return current. As a result, even a small hole in the shield causes enough of a disruption to produce a sensible 
effect that can be measured at either end of the cable. The focus on shielded cable types is therefore crucial, 
as this represents the most viable avenue for successful chafe fault detection in the field. 

An additional advantage of shielded-twisted-pair cable, for instance over coaxial cable, is that it provides 
two tranmission-line modes — common and differential — for signal propagation. Since the differential mode 
is used almost exclusively for the transmission of information or power over the cable, the common mode is 
available and perhaps more suited for signals designed to assess the health of the cable without affecting the 

* Absent sufficient preventive diagnostics and prognostics tools for wiring, regular inspections and built-in redundancy are 
the only currently available protection mechanisms against critical system failures. 


9 



primary purpose of the cable itself. For this reason, the experimental results presented here are focused on 
the common mode, while the theory is derived for both. 

The detailed physics-based model developed here is important for at least two reasons. First, it enables 
one to simulate fault scenarios and noise conditions in order to study and characterize optimal cable and 
interrogation signal design for fault detection. Second, it has produced an extremely accurate and com- 
putationally efficient model suitable for the design of optimal fault detection algorithms, along with the 
quantification of one’s ability to detect faults under a large variety of practical noise conditions, the details 
of which will be the topic of a future publication. 

1.2 Literature survey 

An extensive survey of previous work in this area (including literature in German and Russian) was conducted 
as part of this study, whose pertinent findings will now be summarized. 

The two-wire (or twin-lead) transmission line, consisting of a pair of parallel round conductors, has 
been used extensively in power distribution and telephony applications, and therefore its study garnered 
much attention in the late 19 th and the early 20 th centuries. Various techniques ( e.g ., method of images [4, 
pp. 162-170], conformal mapping [5, pp. 76-78], variational principles [6, pp. 273-279], etc.) may be used to 
analyze the electromagnetic field distribution around this line. (The requisite presence of insulation around 
the wires — typically a thin dielectric ribbon which helps to avoid shorting and to maintain uniform spacing 
between the conductors — changes the physical picture only slightly, and is therefore typically ignored for 
analytical convenience.) 

The introduction of additional wires, or a shield that envelopes the wires, presents significant analytical 
challenges, which consumed considerable effort. The fundamental difficulty here stems from the fact that 
there is no known coordinate system in which all the conductors can be made to coincide with constant- 
coordinate surfaces. The techniques mentioned above were generalized to these problems, but the solutions 
obtained were now approximate [7, 8, 9, 10, 11]. Perhaps a more principled approach to this problem 
(see, e.g., [12, pp. 38-40 and 53]) entails the representation of charge density around the periphery of each 
conductor as a Fourier series with unknown coefficients, using the two-dimensional free-space Green function 
to express the electric field in the space between the conductors as a sum of integrals over the conductor 
surfaces, and finally imposing the boundary condition of vanishing tangential electric field at each conductor 
surface. This procedure leads to an infinite set of linear equations for the unknown Fourier coefficients, 
which may be solved to any desired accuracy. This approach was used by J. Craggs and his followers in 
calculating the capacitance of the shielded-uniform-pair cable [13, 14, 15, 16, 17, 18]. Unfortunately, this 
solution approach does not yield a closed-form expression for the field distribution inside the insulator, 
which one needs when studying the consequences of deviations from the idealized geometry and/or infinite 
conductivity of the wires and/or the shield. 

For long-distance transmission of signals or power, the main concern is attenuation due to energy loss in 
the conductors. In studying this aspect, metals may no longer be assumed perfectly conducting. Indeed, one 
must now solve for the fields everywhere in the cable, and stitch the solutions together at the metal insulator 
interfaces using the appropriate continuity conditions for the fields. The solutions are again in the forms of 
infinite series whose coefficients may be obtained, at least in principle, to arbitrary orders of accuracy. For 
the sake of historical accuracy, it should be noted that J. Carson pioneered this general approach in his study 
of the proximity effect in cables, and in fact suggested in [19] — some 25 years ahead of J. Craggs — that it 
can be used for the calculation of the cable capacitance (see also [20, 21, 22, 23, 24, 25, 26, 27, 28]). Craggs 
and his followers seem to have been unaware of this earlier work. 

As the density of communication lines increased, cross-talk among cables became a paramount issue. It 
was recognized that the simple device of twisting the wires around each other endowed the two- wire line with 
considerable degree of immunity to electromagnetic interference. As one would expect, analytical difficulties 
persist when one moves to this twisted geometry. The geometric complications arising from the twist may 
be removed to a large extent by working in a curvilinear coordinate system that twists with the wires, and 
the resulting held equations may be approached exactly or via perturbation theory [29, 30] (see also [31, 
pp. 96-101]). Alternatively, the Craggs approach outlined above may be generalized to three dimensions 
by employing dyadic Green functions (see, e.g., [32, pp. 23-25], [6, pp. 121 -130], or [33, pp. 544-549]); 
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the resulting formulation is known as the electric-field integral equation in the engineering literature (see, 
e.g., [34, pp. 696-707], [35, pp. 354-356], [6, pp. 139-153], or [33, pp. 617-624]). In any case, the fundamental 
obstacle posed by the non-separable nature of the conductor surfaces remains, and the final solution can 
only be obtained numerically, as in [36]. 

Instead of the continuous Fourier representations of the surface charge distributions adopted by Craggs, 
the wires may also be represented by discrete volume charges within the bulk of the metal, corresponding 
to infinitesimal current filaments in three dimensions. This approximation — a perfectly valid alternative 
paranretrization of the normal electric field at the conductor surfaces — considerably improves the analytical 
tractability of the problem. Such filament models, wherein the wires are treated as conductors of infinitesimal 
cross-section, were preferred, especially in the early days of modest computing power, in analyzing the 
radiation characteristics of twisted-pair cables [37, 38, 39, 40, 41, 42]. The geometries considered therein 
naturally lacked a shield. Mention must also be made of an excellent paper on scattering of plane waves from 
a uniform pair of round conductors [43], which combines filament and wire approaches nicely, and provided 
inspiration for the present work. 

On a different front, in the 1940s and 50s, designs featuring helically wound conductors became popular 
for high-frequency antennas [44, pp. 265-339] as well as for slow- wave devices such as traveling- wave tubes [45, 
pp. 19-48, 229-232], and considerable modeling efforts were devoted to the understanding of their propagation 
characteristics. Unlike in the transmission-line problem, the field distributions had to be known accurately 
here, since they directly determine the radiation pattern of an antenna or the efficiency of interaction with a 
charged-particle stream. Of particular interest were the sheath helix and the tape helix models wherein the 
current-carrying surface is confined to a shell of infinitesimal thickness [46, 47, 48, 49, 50, 51], with ad hoc 
analytical devices employed to account for finite wire radius [52] (see also [45, pp. 34-43] or [53, pp. 306-308]). 
Textbook overviews of these efforts may be found in [54, pp. 40-45, 46-54, 77-82], [55, pp. 476-478], or [6, 
pp. 637-640] . Since the conductor surfaces of this geometry are separable in the cylindrical coordinate system, 
much progress could be made analytically in this case, especially toward obtaining the field distributions 
established around the helical conductor(s), though a fully satisfactory treatment eluded the early workers 
and did not appear until recently [56]. 

In addition to the field-theoretic models of uniform and twisted cable structures described above, it 
is possible to approach the problem from a circuit-theoretic viewpoint when one is interested only in the 
transmission-line modes of the cable. The twisted-pair cable, in particular, received this type of treatment 
extensively, ranging from a first-principles derivation of the equivalent-circuit equations [57] to the experi- 
mental characterization of equivalent-circuit parameters [58], the application of cascaded T-matrix approach 
on piece-wise uniform models whose parameters were calculated analytically [59, 60, 61] or numerically [62], 
finally to the treatment of finite wire radius by a numerical solution of the electric-field integral equation via 
the finite-element method [63]. This basic modeling approach was generalized recently to shielded- twisted- 
pair cable and validated experimentally [64]. 

1.3 Report overview 

The development in this report blends and refines various elements of the analytical approaches mentioned 
above in an effort to obtain a highly accurate yet closed-form model of wave propagation on shielded-pair 
cables. The report consists of three main chapters. In §2, a physics-based model is developed for signal 
propagation along a cable consisting of two wires enclosed by a metallic shield. Specifically, the eigenvalues 
and the corresponding field distributions for the common and the differential transmission-line modes are 
derived. Amplitudes of forward- and backward-propagating modes are then obtained by considering sources 
and loads attached to the cable. The case of uniform wires is studied in some detail, and the results are 
subsequently extended to the case of twisted wires. 

In §3, a careful series of first-principles simulations is undertaken to verify the theory of §2. By solving 
the exact electromagnetic field equations numerically for the same cable geometry, the assumptions made in 
obtaining the approximate analytical model are checked for validity and accuracy. Alternative geometries, 
which cannot be treated analytically, are also simulated to increase confidence in the applicability of the the- 
oretical model to practical cable constructions. Finally, parameters for a preliminary fault model considered 
in §4 are obtained by simulating the effects of a hole in the shield. 
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In §4, a comprehensive framework is developed for extracting cable and fault parameters from experi- 
mental data. At the core of this framework is the S'-parameter formalism for describing a multi-port system. 
Specifically, the theory of §2 is used to derive the S'-parameters of nominal and faulty cable segments. Con- 
siderations of source and load effects and other features of the experimental setup lead finally to a combined 
measurement system model. The estimation of unknown cable and measurement device parameters from 
laboratory data is then couched as a Bayesian inference problem. The power of the overall modeling and 
inference approach is demonstrated on data from both nominal and faulty cables. 

The main conclusions of the report are collected in §5. 
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Chapter 2 

Theory 


2.1 Basic model formulation 

2.1.1 Cable specification 

Consider the cross-section of an ideal shielded-pair cable, oriented along the z axis of a cylindrical coordinate 
system {p,<f>,z), as shown in Fig. 2.1. The shield and the two (identical) wires are labeled “conductor 0,” 
“conductor 1” and “conductor 2,” respectively, the former having (outer) radius r s , and the latter ?' w . The 
space between the conductors is filled entirely with an homogeneous, isotropic insulator of radius ?'j, having 
permittivity e, permeability /x, and conductivity a. The centers of the two wires are placed symmetrically 
on a circle of radius q. Contact between the conductors is avoided via the constraint 2 r w < q + r w < r t . The 
cable is insulated from the ambient medium by a jacket on the outside of the shield (not shown). 


t! 



Figure 2.1: Shielded-pair waveguide cross-section. 
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If this cross-section is translated uniformly in the positive z direction, one obtains a shielded-uniform- 
pair (sup) cable. If, during such translation, one simultaneously rotates the cross-section at a constant rate 
k (rad/m) in the positive <fi direction, a (right-handed) shielded-twisted-pair (stp) cable results, having a 
period p = 2n/n along the z axis. (This may be termed a “polar” twist; see [65] for a “toroidal” twist.) 
Figures and plots in this chapter are generated for a “typical” STP cable whose geometric and constitutive 
parameters are listed in Table 2.1.* 


Table 2.1: Shielded-pair cable parameters. 



Shield 


Wires 

Insulator 

Parameter 

Value 

Parameter 

Value 

Parameter 

Value 

r s 

1.5 mm 

r w 

0.3 mm 

n 

1.35 mm 



V 

30.0 mm 

e 

2.3 eo 



q 

0.675 nun 

M 

1.0 mo 

cr s 

5.998 10 7 S/m 

CT\N 

5.998 10 7 S/m 

a 

0.0 S/m 


2.1.2 Maxwell equations 

With all the conductors initially assumed to be perfect, the electric and the magnetic fields within the 
insulator satisfy the (source-free) Maxwell equations 

V x E = iw/iH, (2-1) 

V x H = (a - iwe)E. (2.2) 

A suitable basis and a coordinate system must be chosen in which to solve these equations. Among the two 
closest candidates, the polar and the bipolar coordinate systems [66, pp. 1175-1215], the former is preferred 
as it affords a more accurate representation of the fields near the shield. (Since the ultimate goal of the 
present work is to develop a physics-based algorithm for detecting faults in the shield of the cable, fields in 
the vicinity of the two conductors need not be modeled with great precision.) It then becomes necessary, as 
argued in §1.2, to approximate the wire cross-sections with a shape that is more amenable to analysis. It is 
judged that a pair of filaments, helically wound on a (fictitious) cylindrical shell of radius r c £ (q -r w ,q + r w ) 
and with period p, comes closest to mimicking, both qualitatively and quantitatively, the fields produced 
by the round wires. (Another possible choice, a pair of similarly wound tapes with “optimally chosen” 
thickness and width, has sharp features that produce large localized fields, which are uncharacteristic of the 
round wires.) As will be demonstrated in §2.2.3, a judicious choice for the filament radius helps retain the 
salient features of the proximity effect, which may be suspected, at first glance, to be thrown out by this 
approximation. Denoting by ij(z) the current flowing on the j th filament, the (surface) current density on 
the shell at p — r c containing the filaments may thus be written 

K(0, z) = — + Z [h(z) 5(<j) - kz) + i 2 (z ) 6(<l> - kz - tt)] , (2.3) 

t c v 1 + {nr c y 

where S(-) denotes the Dirac delta function. This “surrogate” cable model is depicted in Fig. 2.2. 

Taking the curl of (2.1), substituting from (2.2), and dot-multiplying by z, one obtains the Helmholtz 
equation 

V 2 E|| + fc 2 £)| = 0, (2.4) 


*The exact vacuum permittivity and permeability are eo = 8.85418781762 10 -12 F/m and fio = 1.25663706143 10 -6 H/m, 
respectively. For frequencies in the microwave regime and below, the permittivity and the permeability of (nonmagnetic) metals 
are typically taken to have these values. 
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Figure 2.2: A pair of helically wound, current-carrying filaments inside a cylindrical shield (not shown). 


where k = yfiu;/i(c7 — iwe) is the wavenumber in the insulator, and || denotes the axial component. Floquet’s 
theorem dictates that the fields in an helically symmetric structure must have the form (see, e.g., [67, p. 584]) 


W= E 


i(p) e ’ 


i m(4>— kz) giAz 


(2.5) 


Inserting (2.5) into (2.4), the radial dependence of the axial electric field is found to obey the Bessel equation 


den 


1 d 


p dp \ dp 


m 


= 0, 


( 2 . 6 ) 


where 7 m = -^/(A — mn ) 2 — k 2 . This definition of the radial eigenvalue reflects the expectation that the 
phase velocity for the principal modes of interest should be slightly less than that in the unbounded insulator. 
Elimination of the electric field from (2.1) and (2.2) leads to a similar representation for the axial magnetic 
field. 

Next, cross-multiplying (2.1) and (2.2) with z on the left, one obtains 


9Ej_ 

dz 

<9H| 

dz 


= Vj_£’|| — iwpz x Hj_, 

= — (a — i toe) z x Ej_, 


(2.7) 

( 2 . 8 ) 


where _L denotes the transverse component. Substituting Floquet expansions of the form (2.5) for all the 
field components, (2.7) and (2.8) are transformed into 

, , de|| m im . „ 

i(A - mn)e ±m = p — b 0 — e \\m - iw/rz x hj_ m , 

dp p " 


i(A - mn) hj_ m = p - 


dh\\ 


dp 


X 1171 u 

0 — h \ | 

p 


— (a — iwe) z x ej_„ 


Solving for e± m and hi m , one finally obtains 


— 


1 

/ m L 


I ■, dh|| m » im A ( „ de|| m ~ im 

luip ( 0 — p — h \\ m ) -i(A -rat) ^^+0y e II 


h I m. n 


^2 
/ m L 


(a — icje) ( 0 


deii m „ im A dhii m - im 

dp ^ ~p e H m y — 1 A — mK ) + " 


(2.9) 

(2.10) 


giving the transverse field components in terms of the axial ones. 
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2.1.3 Field solutions 

For the “interior” region (0 < p < r c ), the requirement of finiteness at p = 0 leads for (2.6) to the solutions 


e || m(p) — d m J m (7 m /5), (2-11) 

h^ m ( P ) = BZImilmP). (2.12) 


For the “exterior” region (r c < p < r-,), on the other hand, the requirement of vanishing axial and azimuthal 
electric fields, and therefore ey m and dh\\ m /dp , on the shield at p = r\ leads to 


e Tm(p) = A n 

h Tm(p) = K 


ImiPfmP) 


Im ) 


Kmip/mP) 


li'lmP) / , Km\fmP) 

^mV1mr\) 


(2.13) 

(2.14) 


In these solutions, J m (-) and K m (-) denote the modified Bessel functions of the first and the second kinds, 
respectively, while ' signifies differentiation with respect to argument. 

The unknown coefficients in (2.11) through (2.14) are now determined by imposing the continuity condi- 
tions 


p x (E ex - E in ) = 0, (2.15) 

p x (H ex - H in ) = K (2.16) 


across the current shell at p = r c . Toward this end, (2.3) is put into a more useful form by defining the 
helix pitch angle ip = cot -1 (Kx), substituting ij(z) = Ij e lXz to endow the mode currents with the same axial 
dependence as the mode fields they support, and making use of the identity S((p) = l/(27r) oo eim ^ 1° 

acquire the Floquet form, resulting in 


K ( 0 ,*) = 


4> cos ip + z sin ip 
27T7V 


E [h + {-i) m h] 


) iAz 


The application of (2.15) leads immediately to 


A 1 " = A* 


B'l = Bt 


1 - 


1 - 


dm (7 m r\) K m {l mX c) 


BmipfmT i) 4 m (y m r c ) 


K 'mhmr\) I'mipimrc) \ ' 


The application of (2.16), with the use of (2.17), (2.18), and (2.19), then yields 

4 ex Km(lm,r\) T Tn(\ TTIK) COS Ip + 7m sin Ip r r ^ 

A. m = i m ( 7 m r c ) 7ZT77 T77; Bi + (“1) t 2 J , 


B ‘ ; x = 


lm(7m^i) 
K m {'Ym r i) r / 


2irr c (a — iwe) 


r/ / „ X 7^ cos ip 

^ml7m r c) 0 [-^1 H" ( 1) -^ 2 ] 5 


4(7™n) 2 ^ 

where use was made of the Wronskian I' m {x) K m (x) — I m (x) K' m (x) = 1/a;. 


(2.17) 

(2.18) 

(2.19) 

(2.20) 
( 2 . 21 ) 


2.1.4 Eigenvalue equation 

Finally, the heretofore unknown axial eigenvalue A must be determined by imposing the remaining boundary 
condition of vanishing tangential electric field on the two filaments; viz., 

(■ <p cos ip + z sin^/i) • E = 0 (2.22) 
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at p = r c , <t> = nz and <f> = kz + tt. Substituting from (2.9) through (2.12) and then from (2.18) through 
( 2 . 21 ), one obtains from ( 2 . 22 ) the linear system 


where 


E 


1 (-l) m 

1 

1 — 1 

1 

1 

t-H 

1 

1 


S m (u, A) = 0, 


Sm(oj, A) 





KmiPfmT i) 

d-rn (imT ) 


+ O c ) 2 [I'milmTc )] 2 


. r m{lm.r\) 


I'm(lm r c) _ 


The two nontrivial solutions of (2.23) are readily identified; viz., 


K m (imtc ) 
^m(Tm^c) 


(2.23) 


(2.24) 


I 2 = Ii => even (or longitudinal) mode, 
I 2 = —Ii => odd (or balanced) mode. 


The eigenvalue equation, from which the permissible values of A are to be solved for, is therefore given by 


OO 

A(uj,\)= J2 [1 ± (~ l) m ] S m (uj, A) = 0, 

m =— 00 


(2.25) 


where the top (bottom) sign corresponds to the even (odd) mode. Quantities pertaining to the even and the 
odd modes will be indicated throughout by the subscripts e and o, respectively. 


2.2 Uniform cable 


2.2.1 Eigenvalues 

In the limit p — > 00 , the waveguide geometry reduces to that of SUP cable. Setting k = 0, (2.25) becomes 


A(u>, A) = 


w -1 


E 

m =— 00 


[l±(-in [Imhr c )Y 


K m inr\) _ K m (^r c ) 
I m {ir c ) J 


= 0 , 


(2.26) 


f2(cj,A) 


where 7 = \/\ 2 — k 2 . Evidently, A = ±k satisfy the eigenvalue equation for both the even and the odd modes 
of SUP cable. As will be shown in §2.2.2, the mode fields corresponding to these eigenvalues are transverse- 
electromagnetic (tem). These are, of course, none other than the two degenerate transmission-line modes of 
this uniform three-conductor structure. The remaining roots of (2.26), due to the zeros of the factor <?(w, A), 
correspond to mode fields that are transverse-magnetic, transverse-electric, or possibly hybrid in nature. The 
cut-off frequencies of these waveguide modes may be obtained as the zeros of the function f7(w,0), which 
is plotted against linear frequency / = u>/{ 2n) in Fig. 2.3 using the parameter values in Table 2.1 and the 
optimal filament radii to be found in §2.2.3. Specifically, it is observed that the lowest-order even (odd) 
mode is cut off below roughly 67 GHz (103 GHz). Clearly, none of the waveguide modes can propagate over 
the spectral region of typical operation for a cable, and therefore they need not be studied here. 


2.2.2 Mode fields 


With ip = 7t/2, one has from §2.1.3 that 


ey m(p) = I m {lP<)Im{lP>) 

h\\m{p) = 0 ? 


'KrJjTi) 

. An (in) 


K m (rfP>y i 2 i 

Im{lP>) J 27r(cr - i we) 


[1±(-1H, 
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(2.27) 

(2.28) 



(a) ( b ) 




Figure 2.3: Cut-off frequencies of the even (a) and the odd ( b ) waveguide modes of SUP cable. 


where /?< (p>) denotes the lesser (greater) of r c and p , and I is the current flowing on filament 1. In the 
limit 7 — > 0, small-argument approximations for Bessel functions [68, p. 375] reveal that e\\ m also tends to 
zero, leading to Eu = H\\ = 0 everywhere; i.e., the even and the odd modes corresponding to A = ±k are 
both TEM. 

Now, with (2.28) holding identically, the transverse fields (2.9) and (2.10) become 

, , iA dey m ~ i m \ 

£j_m(p) n ( P A "b 0 ^|| m I 7 

7 2 V d P P ) 

, , s a - live ( ~ ± dei| m „ i m \ 

h_L m\P) 9 ( 0 P ^|| m I • 

7 V dp p ) 


Substituting from (2.27) and subsequently into the Floquet series for transverse fields, one obtains 

E ± (r) = -V ± $(r ± )e ±ikz , (2.29) 

Hi(r) = ±i ff ~ M zxEi(r), (2.30) 

where, passing to the limit A — > ±k using small-argument approximations for Bessel functions and making 
use of the expansion ln(l + x 2 — 2a; cos c/i) = —2 X)m=i( a - m / TO ) C0S ( TOl ^)) valid for \x\ < 1, 


^(rj_) = -- 




In 


r \ + ( 7'cp ) 2 — 2rfr c pcos(j) 


rf + ( r c p ) 2 + 2r?r c pcos<j) 
47 t(<j — iwe) [ r i 2 (P 2 + — 2r c p cos <£) J [ r ? (p 2 + r c + 2 r c p cos (j)) 


± In 


(2.31) 


F(r±) 


It should be noted that the upper (lower) signs in (2.29) and (2.30) refer to forward (backward) propagation, 
whereas those in (2.31) refer to the even (odd) mode. Aside from the prefactor, the function <P is, of course, 
nothing but the (electrostatic) potential due to two charges of same or opposite polarity fixed symmetrically 
on a circle of radius r c inside a metallic cylinder of radius ?’i (see, e.g., [12, p. 90]). Noting that / changes 
sign with the direction of propagation, it is observed that the field components obtained above conform to 
the sign convention of the standard normal-mode formalism (see, e.g., [12, p. 390]). 

The potential distributions for the two TEM modes are shown in Fig. 2.4, plotted using the optimal 
filament radii to be found in §2.2.3. The filament current I is chosen (arbitrarily) to produce unit average 
potential on conductor 1 in both cases (see Fig. 2.8). Consequently, the average potential on conductor 2 is 
+1 V (—1 V) for the even (odd) mode, while the exact potential on the shield is zero in both cases. Note 
that the potentials are plotted only over the insulator, since (2.31) is not valid within the wires. 

The azimuthal distributions of charge on the conductor walls reveal the extent to which the proximity 
effect may be expected to play a role in the subsequent determination of equivalent circuit parameters for 
the cable, as detailed in §2.3. On the j th conductor at, say, z = 0, one has the surface charge density 


Zj = j • VxW), 


(2.32) 
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( a ) 


(b) 



Figure 2.4: Potential distributions for the even (a) and the odd (6) modes of SUP cable; the vertical range 
is [0, 1] V in (a) and [—1,1] V in (6), with shades of red (blue) indicating positive (negative) values. 


which may be evaluated analytically by direct substitution from (2.31) into (2.32), preceded, where necessary, 
by a change of coordinates via (2.35) and (2.36). These charge densities are displayed in Fig. 2.5, using the 
same value for / as in Fig. 2.4. It is observed that the charge densities are highly nonuniform around the 
conductor walls. 

(a) ( b ) 




Figure 2.5: Surface charge distributions on the shield (a) and on conductor 1 ( b ) for the even (red) and the 
odd (blue) modes of SUP cable. 
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Substituting (2.31) into (2.29), the corresponding electric field is readily calculated, obtaining 


Ej_(r) 


ik\I\e ±ikz 
2n(a — i we) 



+ 4 > 


rfp — rfr c cos <j> 

p — r c cos 4> 

r f + (r c p) 2 — 2rfr c p cos (f> 

p 2 +rf — 2 r c p cos (j> 

rf p + rfr c cos <j> 

p + r c cos </> 

rf + ( r c p ) 2 + 2rf r c p cos (j) 

p 2 + rf + 2 r c p cos <f> 

rfr c sin <j> 

r c sin (j) 

rf + ( r c p ) 2 — 2rp r c p cos <j> 

p 2 +rf — 2 r c p cos <fi 

rfr c sm(j) 

r c sin 4> 

rf + ( r c p ) 2 + 2rfr c pcos(/> 

p 2 + rf + 2 r c p cos </> 



(2.33) 


The (rectangular) electric-field components of the two inodes are shown in Fig. 2.6. The magnetic field is 
related trivially to the electric field via (2.30), and need not be exhibited explicitly. 


2.2.3 Choice of r c 

The optimal radius r c for the filament pair in Fig. 2.2 will now be determined by imposing the “proper” 
boundary condition 

hj x E(Wj-) = 0, (2.34) 

rather than (2.22), on the round wires, where W j denotes the bounding wall of the j th conductor with 
outward unit normal hj (see Fig. 2.1). In view of (2.29), (2.34) implies that must be constant on W^; the 
value of r c will be chosen to achieve this as closely as possible for each mode. 

Owing to the symmetry of the cable cross-section, it is sufficient to impose (2.34) on only one of the 
conductor surfaces. With (g,p) denoting the local polar coordinates around conductor 1 (see Fig. 2.7), one 
has that 


p 2 = g 2 + q 2 + 2qgcosp, 
p COS (ft = q + QCOSLp. 


Substituting into (2.31) and subsequently expanding into a Fourier series in ip, one obtains 

2 oo r 


= - 


m 


47 r(cr — iu;e) 


In 


-E 


r c g 


± - 


rf — qr c 
r c g 


r c - q 
0 


rf + qr c 


T - 


r c + q 


cos(nup ) 


valid for r w < g < min(g r , ?’i — q), where 


R e 

Ro 


rf - (gr ce ) 2 
(r ce + q)r 2 ’ 

(r co + q)(rf - qr co ) 

(rf + qr co ) 


(2.35) 

(2.36) 


(2.37) 


(2.38) 

(2.39) 


Among possible strategies for choosing r c , setting the amplitude of the first (azimuthal) harmonic of 
( I > (r w , <p) to zero gives arguably the most satisfying result. This strategy leads from (2.37) to a quartic 
equation for r c , whose analytic solution is too unwieldy to write down (see, e.g., [69, p. 23]); approximate 
solutions, accurate to first order in the small quantity (r c — q), are found to be 

= r w( r t ~~ 3q 4 ) + 2q 2 (rf — q 4 ) 

2 q(r? ~q 4 - 4 r^q 2 ) 

= r w( r f + q 2 ) 2 — 2g 2 (rf — q 4 ) 

2 q [r^ (3rp + q 2 ) - rf + q 4 ] ' 
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(a) 


(b) 



Figure 2.6: Electric-field distributions for the even (a and b ) and the odd (c and d) modes of SUP cable; the 
vertical range is [—5, 5] kV/rn in all the plots. 
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Figure 2.7: Optimal filament positions for the even (red) and the odd (blue) modes of SUP cable; the dashed 
vertical lines indicate the centers of the two wires. 


(a) 


(b) 



c p (rad) 


<p (rad) 


Figure 2.8: Variation of the approximate potential (2.37) on the periphery of conductor 1 for the even (a) 
and the odd (6) modes of SUP cable. 


For the parameter values in Table 2.1, the above expressions give r ce = 0.773523 mm and r co = 0.680362 mm, 
while the corresponding exact values found numerically are r ce = 0.760679 mm and r co = 0.680347 mm. 
These latter filament positions are indicated in Fig. 2.7. With the (exact) numerical values of r c inserted 
into (2.37), the variation in the potential on Wi, defined as the ratio of the amplitude of the second harmonic 
to In (R/ g) 2 at g = r w 4 is found to be around 3 % (see Fig. 2.8). Other conceivable strategies — e.g ., choosing 
r c to jointly minimize the first two harmonics, introducing a second pair of filaments, etc . — do not yield an 
appreciable improvement over the simpler strategy adopted here. 

2.2.4 Exact potentials 

Since V • (V x V) = 0 for any vector V, it follows from (2.2) that V • E = 0. Substituting from (2.31), it is 
seen that the potential must satisfy the Laplace equation in the insulator; viz . , 

Vi<2>( rjJ = 0 for r ± G A. (2.40) 

In keeping with the convention established in §2.2.2, (2.40) is to be solved under the boundary conditions 


4>(W 0 ) = 0, 

<W) = +1, 

4*(W2) = il. 

This task may be performed numerically via the relaxation method (see, e.g., [70, pp. 1059-1066]). 
Toward this end, a 2.5-/j,m resolution Cartesian mesh with a total of 825, 559 grid points is constructed over 
the insulator. The analytical approximation (2.31) is used to generate initial guesses, and the grid-point 
potential values are then relaxed using the Gauss-Seidel iteration with Chebyshev acceleration. The direction 

tThe amplitudes of the higher-order harmonics in (2.37) are more than 10 times smaller than that of the second harmonic. 
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of grid traversal is alternated between successive iterations. Updates are computed at each grid point as an 
optimal blend of “+” and “x” averages over the nearest neighbors, as described in Jackson [12, pp. 47 50] . 
Iterations are terminated when the largest update becomes smaller than 10~ 3 ; convergence is achieved in 
under 20 iterations for each mode. Results of this computation are displayed in Fig. 2.9 as slices along the 
coordinate axes, showing excellent agreement, particularly away from the two wires, with the approximate 
analytical potentials (2.31). 


(a) 



X 


( c ) 



(b) 



(d) 

1.0 - 
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0.0 
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Figure 2.9: Comparison along the x and the y axes of the approximate analytical (solid curves) and the 
exact numerical (dashed curves) solutions for the even (a and b) and the odd (c and d) mode potentials of 
sup cable. 


2.3 Equivalent circuit model 

2.3.1 Basic variables 

Although the foregoing model of SUP cable in terms of electromagnetic fields is fundamental and self- 
contained, a “reduced” model involving voltage and current affords greater ease in describing the inter- 
connection of the cable with external circuits (see, e.g., [67, pp. 72-96]), as well as allowing for the inclusion 
of the effects of finite metallic conductivity in a relatively simple way. Furthermore, the twist rate in typical 
STP cables is extremely small, and one may therefore reasonably expect that the circuit-theoretic approach 
developed here can also be used profitably in modeling the STP cable. 

Owing to the absence of axial field components for the transmission-line modes, conductor voltages Vj 
and currents ij may be defined uniquely from (2.1) and (2.2) using Stokes’ theorem in the usual way. For 
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the j th conductor,-!' then, one has 


v j( z ) = ~ E(rj_, z) ■ dl. 

(2.41) 

ij( z ) = j H(rj_, z) ■ dl. 

J Gj 

(2.42) 


where P 3 is an arbitrary contour in the transverse plane that originates from the shield and terminates on the 
surface of the j th conductor, and Gj is an arbitrary closed contour in the transverse plane that circumscribes, 
in the positive sense, the j th conductor (see Fig. 2.1). The shield voltage and current are, by definition, 


v 0 (z) = 0, 

io(z) = ~Y^iA z )- 

j 


One could substitute from (2.29) and (2.30) into (2.41) and (2.42) to obtain the conductor voltages and 
currents directly from the field-theoretic model of §2.2, but it is instructive to develop the circuit-theoretic 
model in its own right. In particular, it will be seen below that the same two transmission-line modes 
found before will emerge from this formulation as well. Contact must, of course, be made eventually with 
the field-theoretic model in order to obtain the line parameters defined in formulating the circuit-theoretic 
model (see, e.g., [71, pp. 64-185]). 


2.3.2 Transmission-line equations 

Telegrapher’s equations 

Differential equations governing the propagation of voltage and current along the line may be derived directly 
from the Maxwell equations. Toward this end, note firstly that, with E\\ = iJy = 0, (2.7) and (2.8) become 


<9E_l 

dz 

c>H_l 

dz 


= —iojfiz x Hj_, 

= — (ct — i we) z x E_l . 


Now, differentiating (2.41) and substituting from (2.43) leads to 


(2.43) 

(2.44) 


dvj_ 

dz 



<9E_l 

dz 


■ dl = iw/i 


I H ± -(dlxz) 

3b 


iujFj(z), 


(2.45) 


where Fj is the (per- unit-length) magnetic flux linked by the contour Pj . This, in turn, is related linearly 
to the conductor currents, and may be written in the form 


F A Z ) = Yl l iP '/'(')• 

j' 


(2.46) 


Inserting (2.46) into (2.45), one obtains 


dv 

- = 


(2.47) 


where v and i denote the (2 x 1) vectors of conductor voltages and currents, respectively. Note from 
that the elements Ijj / of the inductance matrix L may be obtained by computing the flux linked by Pj 
the condition 


lim — — 

Az^O Az 




ij" {z’) dz' = Sj 


(A). 


(2.46) 

under 


tin this section, the indices j and j 1 will be understood to exclude conductor 0. 
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Next, differentiating (2.42) and substituting from (2.44) leads to 


dij_ 

da 


on 

dz 


■ dl = — (a — i we) 


E_l • (dl x z) = - (?- - iw) Qj(z), 


(2.48) 


where Qj is the (per- unit- length) electric charge enclosed by the contour C j. This, in turn, is related linearly 
to the conductor voltages, and may be written in the form 


Qj( z ) = 

j' 


c jj’ v j’ ( 2 )- 


Inserting (2.49) into (2.48), one obtains 


= ^( G - iwC) v, 


where the conductance matrix is given simply by 


G = - C. 


(2.49) 


(2.50) 


(2.51) 


Note from (2.49) that the elements Cjj> of the capacitance matrix C may be obtained by computing the 
charge enclosed by Qj under the condition 


lim — — 

/iz->o Az 


rz-\-Az 


Vj" ( a ') da' = Sj 


(V). 


(2.52) 


For the symmetric cable cross-section depicted in Fig. 2.1, the definition of capacitance implies that 
one has Cn = C22 = c s and C12 = C21 = c m , where c s and c m respectively denote the self and the mutual 
capacitances of the two conductors. Consequently, the capacitance matrix takes the form 



It may be deduced from physical considerations that one must have 


c s > 0, 
Cm fa 0, 
Cs T C m 1 0 . 


(2.53) 


Voltage wave equation 

Differentiating (2.47) and substituting from (2.50) and (2.51), one obtains the voltage wave equation 

d 2 v . fa 

d^ “ (l ~ 

The conductor voltages are related linearly via (2.41) to the electric field, which satisfies 


i <jj I LCv. 


c> 2 E 


-k 2 E, 


as seen from (2.43) and (2.44). Therefore, one must have that 



fc 2 l, 


where I denotes the (2 x 2) identity matrix. Consequently, the inductance matrix may be obtained from the 
capacitance matrix simply via 


L = ytteC -1 . 


(2.54) 
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Since real metals are of large but finite conductivity, the two transmission-line modes of SUP cable 
are in fact guasf-TEM. The presence of (small) axial field components renders the definitions (2.41) and 
(2.42) for Vj and ij dependent on the choices for Tj and Gj, which, strictly speaking, spoils the unique 
correspondence between the circuit- and the field-theoretic variables. Fortunately, rigor may be restored by 
assuming only that the volume current density J remains purely axial within the metals. This is not merely 
a convenient analytical device for alleviating the difficulties associated with the circuit-theoretic model, 
but is indeed necessary for maintaining consistency with the filament model, which can only support axial 
currents. Under this assumption, the magnetic field remains purely transverse, while the transverse electric 
field remains perpendicular to the waveguide walls; thus, the uniqueness of definition (2.41) is preserved 
for arbitrary . Meanwhile, the ambiguity in definition (2.42), caused by the appearance of a small axial 
electric field in the insulator, may be removed by agreeing to let the contour Gj collapse onto the conductor 
wall W j. 

Now, with U|| = 0, (2.50) remains intact, while the term Vi£|| in (2.7) leads on the right-hand side of 
(2.45) to the additional term — [Ay (Pj) — E\\ (-Po)L where Pj and Pq denote the end points of the path Tj on 
the j th conductor and on the shield, respectively. Since E\\ is proportional to the axial current density inside 
the metals, these additional voltage drops will also be related linearly to the conductor currents. Accordingly, 
(2.47) is generalized, for the case of finite metallic conductivity, by the introduction of a series impedance 
matrix Z, resulting in 

jp = — (Z — iwL) i. (2.55) 

Differentiating (2.55), substituting from (2.50), and making use of (2.51) and (2.54), one now obtains the 
modified voltage wave equation 


d^v (j 

ai = (i~ iw ) ( zc -^')v- 


(2.56) 


For the symmetric cable cross-section depicted in Fig. 2.1, the definition of impedance implies that one 
has Zn = 222 = z s and 212 = 221 = 2 m , where z s and 2 m respectively denote the self and the mutual 
impedances of the two conductors. Consequently, the impedance matrix takes the form 


• 2-171 

Physical considerations lead one to expect that 

3?[2 x ] > 0, 

S[2 X ] < 0, 

I ~s | A 1 2m | 



(2.57) 


2.3.3 Calculation of C 

Filament model 

It has been established in §2.3.2 that, for an ideal cable with infinite metallic conductivity, knowledge of C is 
sufficient for obtaining the remaining line parameters via (2.51) and (2.54). Based on the prescription given 
in connection with (2.49), the self capacitance may be computed via 


Cs 


— e 



hi • V_L^(rj_)cU, 


(2.58) 


where use was made of (2.32), and the appropriate potential satisfying (2.52) is 


<P{r±) 


^e(u) 

K(Wi) 


^o(u) 

£o(Wl) 


(2.59) 


There is a minor difficulty in using (2.37) to construct (2.59), since Wi is not an equipotential surface under 
the filament model, and hence ^(Wi) are not constant. This ambiguity may be removed by agreeing to use 
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the average value of the potential on "Wi. With this convention, substitution into (2.58), along with hi = q 
and d! = f Q r w dip, leads readily to 


(2.60) 


ln(i? e /r w ) ln(i? 0 /r w ) J ' 


For the mutual capacitance, one may again integrate on Wi as in (2.58), but this time taking the difference, 
rather than the sum, of the ratios in (2.59) in order to satisfy (2.52). The result is 


c m = 7re 


1 


1 


In (R e /r w ) ln(i? 0 /r w )J ' 


(2.61) 


For the parameter values in Table 2.1, (2.60) and (2.61) give c s = 111.477 pF/m and c m = —17.3822 pF/m. 
These theoretical values compare favorably with the results of numerical integration on the numerical po- 
tential computation described in §2.2.2, which yields c s = 111.475 pF/m and c m = —17.3847 pF/m. 


Comparison with other models 

It is worthwhile comparing the accuracy of the analytical capacitance expressions obtained here with those 
found in the literature in order to cement one’s confidence in the present model. One typically finds the 
conductor-to-conductor capacitance specified for balanced (i.e., odd-mode) excitation, and sometimes also 
the conductor-to-shield capacitance for longitudinal (i.e., even-mode) excitation. For the filament model, 
these capacitances may be obtained by substituting from (2.60) and (2.61) into the definitions 

c o = 2 (^s C m ), 

C e = Cs “t“ C m . 

Other well-known formulas are listed in Table 2.2, along with their numerical evaluations for the parameter 
values in Table 2.1. Since the calculation of capacitance is a one-time operation, the simplicity of formulas is 
not a key consideration, and accuracy is the chief criterion by which different models are to be judged. That 
all the formulas are of similar form and complexity reveals an essential commonality among the different 
analytical approaches, as mentioned in §1.2. It should be noted that the observations made below regarding 
the relative accuracies of the various models apply only to the cable geometry considered here, and should 
not be expected to hold generally. 

The results of numerical computation are given to serve as a benchmark, and the filament model is 
included to facilitate comparison. The textbook formulas for the two-wire and the coaxial lines are provided 
for comparison with the analogous cases of balanced and longitudinal SUP cable excitation, respectively. 
(The additional factor of two in the formula for the coaxial line accounts for the fact that c e defined here is 
the capacitance between a single conductor and the shield.) Not surprisingly, both give poor results. 

The Bell Labs formula, due to S. Mead, is widely quoted (see, e.g., [25], [67, p. 117] or [55, p. 250]) and 
was evidently obtained by solving, via successive approximations, J. Carson’s original formulation for the 
proximity effect [19, 22]. Comparison with the King formula for c 0 , however, shows that the additional term 
in the Bell Labs formula is not appreciable. The Philips Labs formulas are given without a derivation or a 
source reference; their crudeness was noted by C. Miller. § 

The other formula based on a filament model in Table 2.2 is the one due to A. Gent. In order to establish 
equipotential surfaces on the two conductors, he placed the first pair of filaments at the bipolar foci, which 
are on the opposite sides of the conductor centers from the optimal filament positions shown in Fig. 2.7. As 
a result, he found it necessary to add five more pairs — one pair outside the shield, with the remaining pairs 
along the x axis within the conductors — yet the accuracy he achieved is still relatively poor. 

King, on the other hand, placed his filaments at the centers of the conductors, and adjusted the positions 
of the image filaments outside the shield to achieve zero potential on the shield. (More precisely, he started 
from an exact formulation a la Craggs, but then made assumptions on the cable geometry that effectively 
reduces his approach to a filament model.) Consequently, his formulas agree exactly with the present model 
if one takes r c = q in evaluating R x in (2.38) and (2.39). As one would expect in view of Fig. 2.7, the 

-The transcription of the Philips formula for c e in [15] contains an inconsequential error, which is corrected here. 


27 



Table 2.2: Comparison of SUP cable capacitance values obtained with different models. 


Source 

e/c 0 



27 re/ 

C e 


Formula 

Value 


Formula 

Value 

numerical computation 


0.992975 



1.35991 

filament model 

l n (Gco +9 

0.992978 

In 

rt~(qr ce ) 2 

1.35984 


\ r vj rf+qr co J 



(^ce+Q)^w^j 


two-wire line 

cosh ~ 1 (q/r w ) 

1.45057 




coaxial line 



21n(ri/r ce ) 

1.1473 

Bell Labs [23] 


0.993252 





Philips Labs [72] 
Gent [7] 

King [73, pp. 34-39] 
Miller [15] 


cosh 1 


r w( r w+ 4 9 2 )( T i -4g 2 ) 
( 4 g 2 r*j ) 2 

2 / 2 2 \ 
q r \ -(<? ~ r w ) 


In 


in 


r w r?+(q 2 — 

q-\-y/ q 2 —r 2 r t f — qy/ q 2 — r * 2 
r w y*. 2 g 2 — r 2 

( 2q r? -q 2 \ 

\ r w r'f+q 2 ) 


0.960697 

0.994949 

0.993252 

0.992988 


cosh 


-l 


4 -K 2 -^) 2 


4qr w r? 


in 


in 


2qr„ 


( r?-q 4 \ 

\ 2 qr„r'f ) 


rf +g 4 r. +g a 



4 4(« n ) 2 ' 

rTuyr 

2 


1 

4 ' 

, 2(,r.i)2(r.4+ ? 4) 

(rf-9 4 ) 2 


H^) 2 i 

( 1+ !F?) 

i+(^) 2 

1 Mir-,) 4 

[ 4 (r 4 — q4)2 


1.40448 

1.43954 

1.35391 


numerical agreement between the King formula and the present model for c 0 is good, while that for c e is less 
than satisfactory. 

The Miller formulas are based on the rigorous approach of J. Craggs described in §1.2;^ the expressions 
given here were obtained by a 2 x 2 truncation of the infinite determinants in [15]. Comparison with the 
exact value for c 0 shows that the additional term in the Miller formula does provide an improvement over 
the King formula. Despite the proper treatment of the conductor boundaries, however, it is observed that 
this approach is not superior in accuracy to the filament model developed here. Going up to a 3 x 3 
truncation might improve the accuracy, but would certainly lead to considerably more unwieldy formulas, 
and is therefore omitted. 

2.3.4 Calculation of Z 

Skin effect 

In order to isolate the contribution of the skin effect, the coaxial arrangement of a single wire inside a shield 
will be treated first (see, e.g., [75] or [76, pp. 545-554]), and the familiar boundary conditions n x E = 0, 
n x H = K at the surface of a perfect conductor will be supplanted by the “impedance” boundary condition 

E\\=CK\\, (2-62) 

leading to the identification of a suitable surface impedance for each conductor. 

Consider, then, a round metallic wire of radius r w and (large) conductivity er w , carrying a current through 
an homogeneous medium with constitutive parameters e, /i, and a. The return conductor is a cylindrical 
tube with inner radius r\ > r w , outer radius r s , and (large) conductivity a s . The wire and the shield 

llThe transcription of the Craggs formula for c 0 in [7], with which Gent compared his result, contained an error; the Miller 
formula derived here supplies the correction. Although B. King and coworkers later picked up the same error in [74], their 
criticism of Gent’s conclusions remains valid. 
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are situated coaxially, and therefore the currents on both conductors are azimuthally uniform. With the 
conductor currents assumed to be purely axial, as has been discussed in §2.3.1, the Maxwell equations may 
be solved by the introduction of an (E-type) Hertz potential iT = z,E(p)e lXz (see, e.g ., [76, pp. 28-32]). 
Thus, inside the wire, one has 


E\\(P) = (£w-A 2 )f(p)e iA *, (2.63) 

A'P 

h 4>(p) = -(o- w - iwe 0 ) e lAz , (2.64) 

where £ w = y / iw/xo(c r w — iweo) — y/EJJIoa^ is the wave number in the wire, and 'E satisfies 

= < 2 - 65 > 

Since A ~ k <C £ w , the solution of (2.65) is approximately 

&(p) = AJ 0 (£ w p), (2.66) 


where A is an unknown constant, and J m (-) denotes the ordinary Bessel function of the first kind. 

At p = r w , one has that = H^, where K signifies an effective surface current density [12, pp. 352-356]. 
It then follows from (2.62) that the ratio E^/H^ gives the AC surface impedance Cw of the wire. Substituting 
from (2.66) into (2.63) and (2.64), one thus obtains 


Cw 


Cw do (CwCy ) 

CTw -A ( ffj f'v'j ) 


(2.67) 


Using small-argument approximations for J m (-) [68, p. 360], one readily finds that 


lim 


bw 


o 27 rr w irr^ a w ’ 


which is the correct (per-unit-length) DC resistance of the wire. Meanwhile, large-argument approximations 
for Jm(-) [68, p. 364] lead to the high-frequency asymptote 


Cw 1 - i I up o 

Inn = , / . 

u-t oo 27T r w 27rr w V 2cr w 


Inside the shield, one instead has the solution 

E(p) = BJ 0 (^p) + CN 0 ^ s p), 

where £s ~ y/iuipo <r s , B and C are unknown constants, and N m (-) denotes the ordinary Bessel function of 
the second kind. The corresponding fields are obtained as in (2.63) and (2.64), giving 

Effp) ~ Cs [A?do(CsP) + CN 0 (t; s p)} e iAz , 

H*(p) ^ <TsZs[BJi&p)+CN 1 {S s p)]e iXz . 

At p = r*i , one now has that K\\ = —H^, since the outward surface normal reverses its direction. In 
forming the ratio —E\\/H^ to obtain the AC surface impedance ( s of the shield, it is necessary to explicitly 
evaluate the quantity C/B. Rigorously, one must write down solutions for E\ \ and H $ suitable for the 
dielectric media in the interior and the exterior of the shield, and match them with the above solutions 
across the three boundaries at p = r w , p = r\, and p = r s to find all the unknown constants along with the 
appropriate axial eigenvalue A. For sufficiently large cr s , the result is simply C/B ~ — Ti(C s r s )/A r i(^ s r s ), and 
consequently 

> _ Cs Jq{£,sV\) AA(Cs r s) ~~ AA(CsC) Jl{£,sVs) m 

“ a 5 Jr (£ s n) m (Csr s ) - AA (Cs n) Ji (&*) ' [ ’ 
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Figure 2.10: Magnitudes (a) and phases (b) of the surface impedances of a round wire (magenta) and a 
cylindrical tube (orange) in a coaxial arrangement; the dashed lines indicate the corresponding low- and 
high-frequency asymptotes. 


Using small-argument approximations, it may be shown that 


lim 

w-v o 2irr\ 


1 

7T (rf - rf )(7 S ’ 


which is the correct (per-unit-length) DC resistance of the shield. At sufficiently high frequencies, where the 
skin effect is fully developed, large-argument approximations provide the asymptotic value C/B — > i, which 
leads to 

Hm _Cs_ = _ Cs_ gp^fen) 1 - i / q^o 

uj—>oo 2nri a s 2?rr i V 2cr s ’ 

where Hm\-) = J m {-) + iAm(-) denotes the Hankel function of the first kind. 

The per-unit-length surface impedances of the wire and the shield are plotted in Fig. 2.10 using the 
parameter values in Table 2.1. It is observed that these impedances are indeed inductive, as has been 
anticipated in §2.3.2. 


Proximity effect 

When two round wires and a tubular shield are assembled into cable form as in Fig. 2.1, the proximity 
effect enters into the picture, modifying the azimuthal distributions of currents on the conductors from the 
uniform case of the coaxial arrangement considered above. This situation is analyzed next by perturbing 
the boundary conditions (see, e.g ., [66, pp. 1038-1064] or [12, pp. 366-368]) on the field solutions obtained 
in §2.2.2, thus supplying the connection between the surface impedances and £ s derived above and the 
heretofore unknown phenomenological line parameters z s and z m introduced in §2.3.2. 

With the Hertz potential now taken as II = z<F(r x ) e 1 , the fields in the insulator are 

E_l (r) = iAVi^(ri)e iA “', (2.69) 

£||(r) = (k 2 -\ 2 )V(v ± )e iXz , (2.70) 

Hj_(r) = -(a - iwe) z x V x ^(r x ) e iAz , (2.71) 

where 

Vi!f + (fc 2 - A 2 ) if = 0. (2.72) 

In the “unperturbed” case of perfectly conducting metals, to be indicated by the superscript (oo), the field 
expressions (2.69) and (2.71) must, of course, reduce to (2.29) and (2.30). Since A^°°) = k, this requirement 
translates into 

f (oo) (r X ) = £$(r x ). (2.73) 
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The application of Green’s theorem to the two potentials t I / and ! p(°°) gives 


Jj </' ViP (oo) - V\P] do = j 


* ATI _„»<-> A 

dn dn 


d l. 


(Since the surface normals have been drawn inward the insulator cross-section A in Fig. 2.1, the bounding 
contour 6 = W 0 U Wi U W 2 must be traversed in the negative sense around A.) Substituting from (2.72) 
and (2.73) and recalling (2.40), this becomes 


(fc 2 


a = 
a Je 


& d@ 


d'A / dn dn 




d<P 

dn 


d l. 


At the metal -insulator boundaries, the axial electric field and the azimuthal magnetic field are continuous; 
therefore, one may use the field solutions inside the wire and the shield obtained above to show that 




= c 


cr — i we 


for rj_ £ C- 


d'A/dn A 2 - k 2 

Rearranging, one finally obtains 

(A 2 — k 2 ) 2 JJ da — (A 2 — k 2 ) J <P JJ- d l + (a — iwe) J f 


d'P d# 
dn dn 


d/ = 0. 


(2.74) 


Two approximations are now made to facilitate the solution of (2.74). In the first place, the unknown 
potential IF is replaced, to lowest order of approximation, by its unperturbed version a choice that 

may be placed on a firm footing via the variational principle [6, pp. 255-259]. Note secondly that the change 
in A from its unperturbed value k is expected to be inversely proportional to (some power of) the metal 
conductivities, which are very large; therefore, (A 2 — k 2 ) 2 <C (A 2 — k 2 ). At the same time, the ratio of the 
area integral in the first term to the contour integral in the second term of (2.74) is estimated to be of the 
order of the cross-sectional area of the guide, which, for typical transverse cable dimensions, is extremely 
smallJI Consequently, canceling out the prefactor in (2.31) that is common to all the terms, and neglecting 
the first term entirely, (2.74) simplifies to 


(A 2 



iwe) 



(2.75) 


Filament model 


The two contour integrals appearing in (2.75) will now be evaluated using the potential obtained in §2.2.2. 
In keeping with the convention established in evaluating (2.58), it is again assumed that, with the optimal 
choice of r c , the azimuthal variation of <P on W j is negligibly small; it is important, however, to introduce 
this approximation after the normal derivative of <P, where it is involved, has been taken. Noting that <P is 
exactly zero on the shield and exploiting the symmetry of the cable structure, one readily obtains 




dg 


r w d(p 


= — 167T In ( — 

. r w 


(2.76) 


where substitutions were made from (2.37), and use was made of cos(mtp) cos {m' ip) dip = ir 5 mm f. With 
substitutions from (2.37) and (2.31), and the additional use of Y^m= l xm = x /(l ~ x )> valid for |ai| < 1, one 
similarly obtains 


3o = 2£ w 


/•2n 

’ dF(g , p)' 

L 

dg 


{-2tt 


r w dip + Cs 


dF(p, - cf) 
dp 


1 2 


r\ d 4> 


o/- 87r ? W + ( 7 c - q ) 2 , 8 tt 

z Sw o / \o ' ss 

r w r i - (r c - q) z r\ 


(1 ± l) 2 + 4r] 


i p—r j 

2 (rf + r 2 )T (rf ~rg) i 

rf — rd 


II Exact integration with 'P = <d°°) shows that this ratio is less than 10 6 in magnitude for both modes. 


(2.77) 
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where, as before, the upper (lower) sign goes with the even (odd) mode. 

Finally, the variational expression (2.75) is compared with the circuit-theoretic expression (2.82) for A, 
leading to the identification y 2 = — e(Jo/Ji). Using (2.80) for y e and x 0 , along with the expressions (2.60) 
and (2.61) for c s and c m , this yields the following equations for the unknown impedances 2 S and z m : 


Zs 4~ — 


%s — 


In {R e /r w ) % 

2 7T Ji e ’ 

ln( Rp/r w ) Jo 
2tt Ji d ' 


These are solved readily with the aid of the integrals (2.76) and (2.77), obtaining 


Cw 


27T r w 
Cw 


1 + 


(?’ce ~ g) 2 

, - (r-ce - q ) 2 


<7co - 9) 2 


- (r co - g ) 2 


27rri 


r 4 _ 4 

' i 're* 


+ 


2 rfrj 0 

.,,4 _ „4 

1 i ' rn 


(j*ce - qf 


(?’co - g ) 2 


2t rr w [r 2 - (r ce - g) 2 r 2 - (r co - g) 2 J 27rri 


+ 


Cs 


2r?r 2 

1 I CO 

? .4 _ r 4 


(2.78) 

(2.79) 


2.3.5 General solution 


Voltage and current 

Returning now to the voltage wave equation, the solution of (2.56) requires knowledge of the eigensystem of 
the matrix ZC, which is sought in the form 

ZCx = y 2 x. 

For the symmetric and reciprocal cable model under consideration, one has, from (2.53) and (2.57), that 


ZC 


Zs^s T ^mCm Z s C m 4“ 

2sCm 4“ ^mCs Z 5 C s 4“ Z m C m 


The corresponding eigenvalues and (normalized) eigenvectors are then 


Xe — \/ {Zs 4~ Z m ) (Cg 4“ C m ) , Xo >/(^ ) (Cs c m ) , (2.80) 

Xe = ^[l if, Xo = ^ [1 -1] T 


As expected, the eigenmodes of this circuit-theoretic cable model are found to be the same even and odd 
transmission-line modes of the field-theoretic model developed in §2.2. 

The solution of (2.56) may now be written down in the form 


v(z) = (V+ e iA * z + V- e _iAeZ )x e + (V+ e iA ° z + V" e - iA ° z )x 0 , 
where, with x = e, o, the axial eigenvalues are 

A x (w) = ^Jk 2 - (j - iw) Xx- 

From (2.55), the corresponding mode currents are then given by 

i(^) = -(Z-ia;L)- 1 ^ 

= ^ (U+ e iAeZ - V- e- iA ^)x e + i- (V+ e iA ° z - V~ e- iA ° z )x 0 , 

where the characteristic impedances for the two modes are defined as 

z ^ = ( 2 rn - z 2 ){c 2 - c 2 m ) + i2ujfie(z s c s + z m c m ) 4- (cu^e) 2 
iA e (c s 4“ c m )[(^s ^m)(cs c m ) icu/ie] 
z ,, = (4 - z 2 ){c 2 - c 2 m ) 4- i2ujfie(z s c s + z m c m ) + (c Jfie) 2 
iA 0 (c s Cm)[(zs 4~ -2m)(cs T ^m) icu^te] 


(2.81) 


(2.82) 


(2.83) 


(2.84) 

(2.85) 
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Figure 2.11: Real (a) and imaginary (6) parts of the even (red) and the odd (blue) mode eigenvalues for 
lossy SUP cable; the dashed line indicates the lossless case. 

(a) (b) 




Figure 2.12: Phase velocities (a) and attenuation lengths (6) of the even (red) and the odd (blue) modes for 
lossy SUP cable. The phase velocity in the (unbounded) insulator is c = l/y/Jie = 1.97677 10 8 m/s. 


For the parameter values in Table 2.1, and using the expressions (2.60), (2.61), (2.78) and (2.79) de- 
rived above for the equivalent-circuit capacitances and impedances, the real and the imaginary parts of the 
transmission-line mode eigenvalues A x are plotted in Fig. 2.11. Also of interest are the phase velocities and 
the attenuation lengths of the modes, defined respectively as 

uj ui 

K[A] = /3’ 

1 _ 1 
3[A] a’ 

these are plotted in Fig. 2.12. Finally, the characteristic impedances Z x of the two modes are plotted in 
Fig. 2.13; at 1 GHz, (2.84) and (2.85) evaluate to Z e = 53.8786 + i 0.11659 H/m and = 39.3419 + 
i 0.0841395 fl/m. It is noted that, except for the magnitudes of the characteristic impedances, there is 
almost no perceptible difference in these quantities between the two modes. 

Mode amplitudes 

Finally, the heretofore unknown mode amplitudes V~ must be determined from a specification of the source 
and the load configurations at the two ends of the cable, whose length will be denoted £. Referring to Fig. 2.14, 
adapted from [73, p. 203], it is assumed that two voltage sources vsi and vs 2 with internal impedances Z$ i 
and Z $2 are connected between conductor 1 and the shield and conductor 2 and the shield, respectively; 
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Figure 2.13: Magnitudes (a) and phases (6) of the even (red) and the odd (blue) mode characteristic 
impedances for lossy SUP cable. 


in the most general configuration, the connection to the shield is made via a common impedance Z$o- At 
the load end, the line is terminated by two impedances Z li and Z \_2 that are connected, on one hand, to 
conductors 1 and 2, respectively, and on the other hand to the shield via a common impedance Z lo- That 
the termination of a transmission line can be modeled in this relatively simple manner may be justified by 
careful analysis, as shown in King [73, pp. 58-73]. 

Writing Kirchhoff’s voltage law around the two separate circuits at both the source and the load ends, 
one obtains the following four equations: 

ui(0) = vsi - Zsizi(0) + ^so*o(0), 

^ 2 ( 0 ) = VS2 — -^ 52 * 2 ( 0 ) + -£so*o(0), 
v\(t) = Z\_iii(€) — Zi 0 i 0 (£), 

V2{£) = Z\_ 2 l2{£) — Z\_ 0 i 0 {t). 


Remembering that io = —i\ — * 2 , and substituting from (2.81) and (2.83), these equations lead to the system 


1 + Asi 1 — Asi 1 + i?si 1 — B%\ 


V 


vsi 

1 + As 2 1 — As 2 — (1 + -0S2) — (1 — i?S 2 ) 


v e - 

= V2 

VS2 

(l-A L i)e iA ^ (1 + A L i) e _iAe ^ (1 -B L1 )e iX ° e (1 + B L1 ) e~ iX ° e 


K + 


0 

(1 -A L2 )e iA ^ (l + A L2 )e- iA ^ -(l-R L2 )e iA ^ -(1 + B L2 )e~ iX ° e 


v 0 ~ 


0 


where, with X = S, L and j = 1, 2, 


( 2 . 86 ) 


Axj 

Bxj 


Zxj + 2Zxo 
Zxj_ 

Z 0 ' 


z = 0 z = i 




Figure 2.14: A transmission-line circuit featuring STP cable. 
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The general solution of (2.86) is readily found but rather cumbersome to write down; the result simplifies 
considerably for the practically interesting case of “balanced” source and load configurations, for which 
Zsi = Zs 2 = Zs and Z L i = Z\_ 2 = Z \_ , yielding 


V 0 ~ 


Z e 1 

Z s + 2Z S0 + Z e l- r Se r Le e i2 ^ 

z 0 1 


Zs + z a i — r So r Lo e 


i2A 0 £ 


1 

Tlo 


1 


r Le 


VS1 + VS2 

V2 


Vsi — VS2 


where the reflection coefficients are defined as 


/ Xe = 


Tko = 


Zx + 2Z\o — Z e 
Z x + 2Zxo + Z e 
Zx — Z 0 
Zx + Z D 


(2.87) 

(2.88) 


Concept of operation 

Note from (2.87) and (2.88) that the even-mode and the odd-mode amplitudes are proportional to the sum 
and the difference, respectively, of the two source signals, a separation that occurs only in the special case 
of balanced source and load terminations considered above. This observation motivates the definition of the 
common - mode and the differential - mode voltages and currents (see, e.g., [77] or [71, p. 37]) 

v c = \ fl + v 2 ), V d =V 1 -V 2 , 

i c =ii+ *2, id = 2 (*i — ii)- 


These definitions are especially useful in analyzing experiments where a single source is used to excite the 
cable (see §4.2). Introducing the new voltage and current vectors 


1 l 

2 2 
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1 1 

1 _1 
2 2 


and substituting from (2.81) and (2.83), one finds that 
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where the new characteristic impedances and basis vectors are 

Z c = \ Z e , Z d — 2Z 0 , 

y c = [i of, y d = [o if, 

and the new mode amplitudes are given by 
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A possible concept of operation is to send a “message” signal v m for communication, while simultaneously 
sending a “test” signal v t down the line for monitoring the health of the cable. Although not discernible 
from Fig. 2.12, the odd mode is the slightly faster and less attenuated of the two transmission-line inodes. 
Note also from Fig. 2.5 that, for the even mode, the current flow on the shield is slightly less nonuniform 
over the azimuth, thus making it more sensitive to shield imperfections such as chafing. For these reasons, 
it seems sensible to assign the message signal to the differential mode, while the test signal is transmitted 
via the common mode. Thus, in the balanced case, the source signals may be taken as 

vsi = v t + \v m , 

VS2 = V t ~ §U m . 


2.4 Twisted cable 


2.4.1 Eigenvalues 

For finite p, the waveguide geometry is that of STP cable. In a typical cable, the twisting of the wires is very 
gradual in comparison with the transverse dimension ( i.e ., p ri), and therefore the eigenvalues and the 
corresponding fields for the modes of STP cable may be expected, at least over the spectral region of interest, 
to behave substantially like those of SUP cable discussed in §2.2. Specifically, for the transmission-line modes 
of STP cable, one expects to find A > k. Furthermore, for the parameter values in Table 2.1, one has that 
k k up to 1 GHz. Consequently, one may adopt the approximation r ) m ~ \m\n in (2.24). Noting that the 
modified Bessel functions and their derivatives obey the reflection relation W- m (x) = W m (x), and making 
use of the corresponding large-order approximations [68, p. 378], (2.25) thus becomes 


A(u>, A) ~ 


k J \/l + (Kr c y 


- \A + ( KT C ) 2 


^ 1 {e- 2 -[H(n)-H(r c) ] — f} = 0, 


(2.89) 


where m runs over positive even (odd) values for the even (odd) mode, and 


“ (p) = V 1 + Op ) 2 + In 


Kp 

1 + yjl + (up) 2 


It should perhaps be noted that, for the even mode, the to = 0 term in (2.25) must be treated separately 
using small- argument approximations for Bessel functions. The contribution of this term becomes ultimately 
immaterial, however, since the sum in (2.89) is technically divergent. The eigenvalue is therefore determined 
by the factor outside the sum, giving simply 

A = ±k \J 1 + (nr c ) 2 . (2.90) 

Remembering that + (kx c ) 2 = 1/ sin '0, it is observed that, at least to this order of approximation, 
the multiplicative “correction” of (2.90) to the eigenvalues of the transmission-line modes of SUP cable is 
purely geometric (i.e., independent of frequency), and essentially accounts for the elongated path the fields 
must travel along the twisted wires in STP cable. It is also intuitively satisfying to note that the correction 
depends on n through an even power, and is therefore insensitive to the sense of the twist (i.e., the sign of 
k), as one would expect on physical grounds. Using the parameter values in Table 2.1, the radical in (2.90) 
evaluates to 1.01261 (1.0101) for the even (odd) mode, revealing that the correction is indeed slight. Note 
that the most prominent feature of periodic structures, namely the forbidden spectral regions known as “stop 
bands” wherein axial wave propagation ceases, is not observed. As was anticipated by early workers [49], 
this feature is eliminated here by the presence of the shield, which traps the radiation that would flow out 
radially in these frequency bands if the cable were an “open” structure. 


2.4.2 Mode fields 

The exact field solutions given in §2.1 may now be evaluated with (2.90) to the same order of accuracy as 
(2.89), making use of both small-argument and large-order approximations for the Bessel functions. Note 
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that it is important to properly account for the m = 0 term here. This procedure leads to approximate 
analytical expressions for the mode fields that are, unfortunately, rather unwieldy to include here. It is 
therefore merely noted that these approximate solutions exhibit a slight twist compared with the field plots 
for SUP cable given in §2.2.2, but are otherwise of the same general character. This result, corroborated by 
exact numerical solutions (see §3.5), was to be expected given the very small twist rate of the wires typical 
for STP cables. 

Finally, it seems reasonable to expect that the equivalent-circuit model for SUP cable may be useful for 
STP cable as well. If the development of §2.3 is followed with the new eigenvalue (2.90) and without regard 
to the small axial field components that appear as a result of twisting, it is found that the per-unit-lengtlr 
line parameters C, G, L, and Z are all increased by the factor 1/ sin^, which is consistent with the geometric 
argument presented above. Consequently, the eigenvalues for the transmission-line modes of lossy STP cable 
may be obtained simply by multiplying (2.82) by the same factor, leading to a small increase in attenuation 
along with the small reduction in phase velocity indicated by (2.90). (The characteristic impedances Z x 
remain the same, as they involve ratios of the per-unit-length line impedances to the eigenvalues.) 
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Chapter 3 


Simulation 


3.1 Overview 

Electromagnetic wave propagation simulations are presented as performed with COMSOL Multiphysics sim- 
ulation software (version 4.2) with the RF and the AC/dc modules. The purpose of the simulations is to 
characterize the cables as a function of frequency in both unfaulted and faulted conditions. This character- 
ization can then be compared with the theoretical derivations as part of a validation process. The results 
can also be used as part of a trade-space assessment for diagnostic tool development for STP cable. 

3.2 Geometric fidelity 

Different representations of the geometry of STP cable were considered. These representations were examined 
in order to determine the greatest level of agreement between the actual cable and the simulator while allowing 
for the computations (or more specifically the memory usage) to be practical. All simulations were performed 
on a dual quad core Intel Xeon W5590 @ 3.33 GHz computer equipped with 192 GB of RAM. 

There are essentially two different kinds of STP cable: impedance controlled and non-impedance con- 
trolled. The difference between these two types of cables is the uniformity of the spacing between the inner 
conductors and the shield. The non-impedance controlled cable can have flat faces and the impedance con- 
trolled cable has been manufactured to be round. From a modeling perspective, we examined four variations 
as shown in Fig. 3.1. The first is the most practical from an analytical perspective, it has a perfectly round 
shield filled with dielectric and two twisting conductors. In practice, the cable actually has air gaps between 
the dielectrics on the wires and the shield. The early days of analytical wave propagation analysis also 
included approximating the inner conductors with infinitely thin metallic tapes. We simulated all of these 
including the non-impedance controlled (flat) configuration. 

The theoretical model had dimensions corresponding to a shield radius of 1.5 mm, a dielectric radius of 



Figure 3.1: Cross-sections of simulated cable geometries. On the left is the geometry corresponding to the 
theoretical model, next has the inclusion of air, and to the right are the typical tape approximation and the 
flat non-impedance controlled cable. 
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Figure 3.2: Meshes for STP cable and hole fault. Top left shows the hole in the shield. Top right shows 
the meshing detail for the hole and the density near the edges. Bottom left shows the mesh, including the 
boundary layer mesh, for the whole 3D structure. Bottom right shows the boundary layer meshes for part 
of the shield and one conductor. 


1.35 mm, and core conductors of 0.3 mm. One 360-degree rotational twist occurred within a 30 mm length. 
The impedance controlled model, which included air, had the dielectric radius around each conductor as 
0.675 mm. All other parameters were equal to the theoretical model analyzed in §2. The electrical parameters 
were consistent across models. The relative permittivity of the insulation was 2.3, and the conductivity of 
the copper was 5.998 10 7 S/m. 


3.3 Mesh 

Most of the 3D simulations presented will rely upon the simplest geometric model. The mesh for this model 
is presented in Figure 3.2. The mesh consisted of four different levels of spatial resolution in order to resolve 
the fields of interest. The perfectly matched layer absorber at the end of the cable is not shown. The 
unfaulted shielded twisted pair had a maximum swept mesh spacing of 1.09 10” 4 m. For the cases with a 
hole fault, the square holes were modeled with a high-density mesh at the very edges for 2 /mi, the rest 
of the volume of the hole was at 40 /mi. A typical 60 mm length of cable without faults used 1.3 million 
elements representing 18.1 million degrees of freedom. The analyses with faults used much shorter lengths 
of cable (typically 5 mm) in order to have enough memory for the iterative solver. These numbers do not 
include when boundary layer meshes were added to facilitate modelling the submicron skin effects for the 
attenuation calculations. The boundary layer mesh shown in Fig. 3.2 (and described later in more detail) 
resulted in using 892,000 elements so that only a 30 mm length of cable for the attenuation analysis was 
able to fit within memory using an iterative solver. 

For the 2D analyses the capacitances were solved for by using a direct MUMPS solver. The 3D eigenfre- 
quency analyses were also performed with a direct MUMPS solver. The 3D field solutions for the cables with 
a shield hole fault were calculated using a stationary iterative solver in the frequency domain based upon 
GMRES with a geometric multigrid. The iterative solver was necessary due to the large memory requirements, 
which were in excess of the machine’s capability for a direct solution. 
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3.4 2D analysis 

The first analysis consists of a 2D (Laplace) solution of the fields in order to determine the capacitance. The 
calculated capacitances are shown in Table 3.1 and the resulting fields are shown in Fig. 3.3. 

The capacitance of the theoretical model and the impedance-controlled model (with air) are within 10 % 
of each other and so it seems that the theoretical approximation is close. The tape model is further away 
from these capacitances and has magnitudes of the electric fields that are significantly peaked due to the 
thin edges. The flat (non-impedance controlled) cable is significantly different as expected. The magnitude 
of the fields is important when we go on to consider the magnitude of the fields in the hole faults. If the 
fields are significantly stronger ( e.g ., tape) than that expected analytically, the expected reflections will be 
significantly in error. The tape model and the flat non-impedance controlled model will not be discussed 
further in our presentation of the simulation results. 


Table 3.1: Capacitances from 2D analysis. 


Model 

Self Capacitance (pF) 

Mutual Capacitance (pF) 

theoretical stp cable 

111.616 

-17.511 

STP cable with foam & air 

95.7642 

-17.8204 

flat stp cable 

257.0246 

-25.1478 

tape 

89.0881 

-5.4137 



Figure 3.3: Pairs of real E x (left) and real E y (right) fields for the even mode. Note that red represents 5000 
V/m and blue represents —5000 V/m. 


3.5 3D analysis 

A full 3D segment representing more than one complete twist was analyzed in COMSOL to determine the 
dominant modes of propagation. An eigenfrequency analysis was conducted on 180 mm of coaxial cable 
(configured as RG-58), on the theoretical STP cable with lengths ranging from 90 mm to 360 mm, 90 mm 
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of STP cable with foam and air (shorter due to increased mesh elements), and 180 mm of SUP cable. The 
computed eigenvalues at the resulting frequencies are presented in Tables 3.2 and 3.3 for the two primary 
modes (even and odd) of the three-conductor cable. The transverse electric fields for these results are 
presented in Fig. 3.4. Table 3.3 modeled the high-frequency skin effect, and so it also includes attenuation. 
It was also confirmed that these computed eigenvalues agree very well with the theoretical model of §2 up 
to — and possibly beyond— 10 GHz. 


Table 3.2: Eigenvalues from 3D analysis. 


Model 

Frequency (GHz) 

Odd-mode 5J[A] (1/nr) 

Even-mode 5ft[A] (1/m) 

coaxial cable 

0.552 

17.36 


coaxial cable 

1.10 

34.56 


STP cable with air 

5.98 

19.01 

19.89 

STP cable with air 

1.79 

57.02 

59.67 

sup cable 

5.49 

17.45 

17.45 

sup cable 

1.10 

34.80 

34.79 


Table 3.3: STP cable eigenvalues from 3D analysis for theoretical model. 


Frequency (GHz) 

Odd-mode A (1/m) 

Even-mode A (1/m) 

0.54 

17.1452 + i0.0522 

17.2476 + i 0.0515 

1.62 

51.4973 + i 0.0809 

51.8045 + i 0.0813 

2.7 

85.8547 + i 0.0958 

86.3679 + i 0.0977 

3.8 

120.2127 + iO. 1056 

120.9327 + i0. 1091 

4.9 

154.5704 + i0.113 

155.498 + i0.118 

5.9 

188.9276 + i0.119 

190.0636 + i 0.1254 

7.1 

223.284 + i 0.124 

224.6293 + i0. 1318 

8.1 

257.6396 + i 0.1284 

259.1952 + i0. 1375 

9.2 

291.9943 + i 0.1324 

293.7612 + i 0.1428 


The transverse electric field resulting from the mode analysis was compared with the solution of the 
Laplace equation in order to assess whether it was representative of the even mode or the odd mode. A 
careful examination of Fig. 3.4 reveals the asymmetrical effects of the twists on the fields in the first pair of 
images as opposed to the last pair that are untwisted. The subtlety of the change in fields between twisted 
and uniform conductors is also revealed by the similarity in eigenvalues between the twisted and untwisted 
cases. The inclusion of air in the impedance controlled model produces a minor change in the eigenvalues. 
(Note that the frequencies reported in the two tables are different due to the inherent differences in velocity 
of propagation and lengths of samples simulated.) 

The calculation of the attenuation in STP cable required significantly more attention to the meshing of 
the skin of the conductors and the shield. Within the mesh generator in COMSOL is a feature to create 
boundary layer meshes. This feature was used to generate a boundary layer mesh specified to have 15 layers 
with a stretching factor of 1.4 and a starting mesh size of 10~ 7 m. This resulted in boundary mesh layers 
at approximately 0.1, 0.24, 0.44, 0.71, 1.09, 1.63, 2.39, 3.44, 4.92, 6.98, 9.87, 13.92, 19.59, 27.53, 38.64 fj.ni. 
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Figure 3.4: From left to right, the first pair of plots depicts the magnitude of the transverse electric field 
for the first two eigenmodes of the theoretical STP cable; the second pair is the impedance controlled cable 
which has air between the dielectrics; the third pair is the SUP cable. All plots are scaled 0 to 50,000 V/m. 



Figure 3.5: Electric field penetration into the metals via the boundary layer mesh of 3.44 to 4.92 /im at 
300 MHz. 


Obviously, the conductors are thicker than 38 //m, so the rest of the volume of the conductors and the shield 
are represented with tetrahedral elements. 

A rough estimate of skin depth ^/2/(w/nr) ranges from 3.75 /xm at 300 MHz, to 2.06 /xrn at 1.0 GHz and to 
1.02 /im at 4 GHz. The dynamic range of the electric field penetrating into the metal was roughly ±10 -3 V /m. 
Fig. 3.5 depicts the electric field within the metal at 300 MHz (the maximum depth of penetration). The 
simulation resulted in depths in the range of 3.44 to 4.92 /tin for 300 MHz, 2.39 to 3.44 /im for 1 GHz, and 
1.09 to 1.63 /iin for 4 GHz. The simulation tended to allow greater penetration because the mesh could 
not be made fine enough to truly represent the skin effect. A coarser mesh forces the fields to exist deeper 
into the metals than usually occurs, but due to computer memory limitations, the mesh density was deemed 
adequate given the penetration ranges. 

The eigenanalysis with this boundary layer mesh resulted in an attenuation of 0.052/m at 540 MHz 
and 0.081/m at 1.63 GHz. A frequency-domain simulation using this mesh was then performed with a 
tem source (x.e., a coaxial port) to obtain the scattering parameters for the common mode. The throughput 
for a £ = 30 mm length of STP cable, given by — In |,SYi |/^, is plotted in Figure 3.6 along with the linear 
loss from the eigenanalysis. The key point to examine in this figure is where S 21 dips to its lowest where 
it approaches the linear attenuation, which matches closely with the attenuation line drawn using the lossy 
part of the eigenvalues. (The oscillatory behavior of S 21 is due to the mismatched terminations at the two 
ends of the cable.) 

The theoretical analytic model compares well with the simulation results as shown by examining the 
transverse fields of the analytical model for SUP cable with the transverse fields of the 2D simulation and 
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Figure 3.6: Attenuation per meter in theoretical STP cable. Shown for a 30 mm length of cable are the 
through scattering parameter — In \S 21 \I and the linear loss line derived from the eigen-analysis. 


the frequency based simulations for the full 3D simulation model. Figs. 3.7 and 3.8 show the profiles of the 
real part of E x when following along the x-axis, and the real part of E y when following along parallel to 
the y - axis through conductor 1. The E x profiles show three field regions from left to right: between shield 
and left conductor, between conductors, and finally between the right conductor and the shield. E x matches 
closely with the analytical model (6 % error) when near the shield and shows the greatest error when closest 
to a conductor (19% error). The most important region for accuracy of detecting a shield hole is the field 
between the conductor and the shield. The E y profiles show two field regions: the fields between the bottom 
of the right conductor and the shield, and the fields between the top of the right conductor and the shield. 
Again, the agreement is close between theory and 2D simulation (1 % error) when closest to the shield, and 
the error is greatest near the conductor (7% error). The error closest to the conductor in both cases is due 
to the filament approximation to the round conductor. The fields at 1 GHz naturally diverge away from the 
lossless 2D solution, and are shown to emphasize that the loss with respect to frequency dominates most 
other effects. 


3.6 Chafe faults 

Since the angular positions of the conductors within STP cable change as a function of length, the field 
intensity rotates with distance along the cable. Thus, when determining fault characteristics, we have to 
consider not only distance to fault and the size of the fault but also the angular position of the fault, which 
can cause substantial variations in reflection depending upon whether it is close to a conductor or far away. 

This key point is illustrated in Fig. 3.9, which plots E y for three discrete hole angles. The faulted 
specimen was 5 mm in length with the fault centered at 2.5 mm. The angle of rotation of the conductors 
was 12 degrees per millimeter. Thus, the center of a hole located on the x-axis at 2.5 mm is at —30 degrees 
with respect to the axis running through the conductors. The other two angles were at —75 and 15 degrees. 
The wave propagation was initiated by specifying E x and E y in the dielectric at one end of the cable. The 
initial values of E x and E y were taken from the mode analysis solution for the primary eigenvalue of even 
mode propagation. 

The fields escaping through the hole are determined mainly by E z in the hole. Fig. 3.10 displays the 
E z field along the centerline of the hole running on the 2 -axis (direction of propagation) for the three hole 
positions. Note that the peak fields at 15 degrees are more than two times in magnitude compared with 
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Figure 3.7: Profile of real E x parallel to x-axis for analytical theory, 2D simulation, 3D simulation at 300 
MHz, and 3D simulation at 1 GHz. 


those at —75 degrees. These peak fields, when adjusted by the area of the fault, are proportional to the 
expected magnitude of reflection from these faults. This implies that it may be possible to have a small 
fault near a conductor be nearly equal in reflection magnitude to a larger size hole fault that is farther away 
from the conductors. These multiple possibilities would then require that any fault parameterization search 
routine be able to cope with such an ambiguity. 

In order to quickly take advantage of the changes with respect to fault angle, a capacitance look-up table 
was created by solving the Laplace equation for the impedance-controlled cable (including air). Note that, 
since this is a 2D calculation, the effects of twist are not captured. Table 3.4 lists the resulting capacitances 
for both conductors with different sized faults at discrete angles. The fault size of 3 mm represents half of 
the shield without intrusion into the dielectric. We could have chosen to also abrade/remove the dielectric, 
but this would have had a greater impact on the results. Thus, by not removing the dielectric, it is possible 
to examine the least detectability inherent within the problem. 


Table 3.4: Capacitances (in pF) with faults for the even mode of impedance controlled STP cable with air 
but no foam. 


Hole size (mm) 

Angle = 0 deg. 

Angle = 45 deg. 

Angle = 90 deg. 

0.0 

71.476, 71.476 



0.5 

71.464, 66.000 

71.436, 70.894 

71.383, 71.383 

1.0 

71.422, 55.682 

71.310, 68.863 

71.066, 71.066 

2.0 

71.181, 37.979 

70.729, 55.760 

69.144, 69.144 

2.5 

70.847, 31.002 

70.062, 42.379 

66.001, 66.001 
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Figure 3.8: Profile of real E v running parallel to y - axis through center of right conductor for analytical 
theory, 2D simulation, 3D simulation at 300 MHz, and 3D simulation at 1 GHz. 
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Figure 3.9: E y field strength at 2 GHz with hole at —30, —75 and 15 degrees from axis formed by conductors. 



z (mm) 

Figure 3.10: Real part of E z along the centerline of hole for three different angular hole positions. 
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Chapter 4 


Experiment 


4.1 Overview 

Before the theoretical model can be compared with experiment, the effects of the measurement setup itself 
must be incorporated to produce a complete model for the entire system that includes instrument impedance 
mismatch, measurement delays, and a small degree of allowable miscalibration, all of which significantly affect 
the measured signal by a much greater margin than the chafe faults we aim to detect. This overall system 
model will now be developed within the ^-parameter framework detailed in [78]. Optimal estimation of the 
unknown or uncertain cable and measurement setup parameters from laboratory data will then be couched 
as a Bayesian inference problem, whose solution will be demonstrated by several examples. 


4.2 System model 

4.2.1 Common-mode cable model 

As discussed in §2.3.5, of the two transmission-line modes of STP cable, the common mode is more suitable 
for the diagnostic assessment of cable health. The pertinent model parameters for common-mode excitation 
of STP cable are the twist-corrected axial eigenvalue 


A c = 


sin ip 


= y [1 + (Kr C e) 2 ] iuj[i(<7 - iwe) - - iw) (z s + z m )(c s + c m ) 


and the characteristic impedance 


£e _ (~m - z s)( c s - C rn) + i2uj/u,e(z s c s + z m c m ) + (ca^e) 2 
2 i2A e (c s T c m )[(zs ^m)(Ts ^m) icc/re] 


(4.1) 


(4.2) 


The line parameters c s , c m , z s , and z m have been derived in §2.3.3 and §2.3.4. In the inference examples 
described in §4.4.1, the geometric parameters of the cable cross-section will be assumed known, and the 
material parameters e, er, <r w , and a s will be estimated from laboratory measurements. In order to explicitly 
exhibit the dependence of the model on the inferred parameters as well as on frequency, the expressions for 
the line parameters are rewritten as 


c s = r s e, 


Cm — 








(4.3) 

(4.4) 

(4.5) 

(4.6) 
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where 1 


K [ln(i? e /r w ) + ln(i? 0 /r w )J ’ 

r 1 i 

ln(i? e /r w ) ln(i? 0 /r w ). 

(w) = [do , ( r ce - g ) 2 ( r co - q ) 2 

5 2tt r w V 2 [ r 2 -(r ce -g) 2 r 2 - (r co - g) 2 J ’ 

(s) = [do ( r t + r i 2r i 2 r c 2 p A 

27r?’i V 2 V 7 f^ r ce r t~ r to)' 

(w) = 1 - i f^o [ (r-ce ~ g ) 2 (r C Q - g ) 2 

??m 2 tt r w V 2 |_r 2 — (r ce — q) 2 r 2 - (r co - q) 2 \ ' 

(s) = 1^2 [*L f r l + _ 2r i 2r co A 

m 27r?’i \ 2 \rf -r£ e rf - ' 


It should be emphasized that, for a given cable geometry, these constants need be evaluated only once. 

For a nominal cable segment of length Az and a very small twist rate, the ^-parameter matrix takes the 
form 


S (n) 


o si n i } 

4i } 0 


where 


c( n ) _ „iA c .4z 
°21 — e 


(4.7) 


4.2.2 A simple chafe fault model 

Deferring a detailed electromagnetic model of scattering from a hole to a later publication, the approach 
taken in [78] for coaxial cable will be followed towards developing a simple fault model. Consider, then, a 
small chafe located at a distance df on a STP cable of length t, having width wf = r\A<f> f, length gf, and 
angular position <fif, as depicted in Figure 4.1. It is assumed firstly that common-mode propagation across 
the chafed cable segment may be modeled by a modified axial eigenvalue Af and a modified characteristic 
impedance Zf , as indicated in the figure, and that the inevitable coupling of energy into the differential 
mode will be negligibly small. It is further assumed that these fault model parameters may be obtained 
from the general expressions (4.1) and (4.2) by recomputing the line parameters therein for the chafed cable 



Figure 4.1: (Left) A constant- impedance model for a chafed cable segment. (Right) Cross-section of a cable 
with a chafed shield. 

*Note that high-frequency approximations are being used here for the line impedances, in anticipation of the frequency 
range of laboratory measurements reported subsequently. 
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cross-section. Specifically, since the (electrical) length of the chafe is very small, the effects of loss due to 
finite conductivity may be neglected by setting a = z s = z m = 0. Replacing the (per-unit- length) even-mode 
capacitance c s + c m with the (per-unit-length) common-mode fault capacitance Cf = Tf€. one thus obtains the 
expressions 


Af = u>vVe[ 1 + (nr ce ) 2 ], (4.8) 



Unfortunately, Tf appearing above cannot be obtained in closed form, and a look-up table must instead 
be constructed off-line. Toward this end, we choose a set of angular fault positions within the range (0, 7t/2) 
and a set of fault widths within the range (0,7rrj). For each possible combination of <j > f and Wf, the 2D 
Laplace equation (Vj_# = 0) is solved under even-mode boundary conditions ( i.e ., both inner conductors 
set to unit potential while the shield is held at zero potential). Since the cable cross-section is now “open,” 
the computational domain must be bounded by an outer surface held at zero potential that is sufficiently 
removed from the shield in order not to disturb the fields leaking out through the hole. The resulting 
potential is then integrated numerically to obtain 


Wf) 


1 

2 


® hi • V_l^(i\l) d l + 
J Wi 


® n 2 • Vj_<£(rj_) dZ , 
J w 2 


(4.10) 


where the dependence of Tf on the underlying fault parameters 4 > f and Wf has been indicated explicitly, and 
the average of separate integrals around each conductor is taken in an effort to compensate somewhat for 
the asymmetry of the cable cross-section due to the presence of the holed 

Since the chafe model is symmetric with respect to the cable axis, the S matrix for the fault with respect 
to the common-mode impedance Z c takes the form 


S (f > = 


q(f) q( 0 
J 11 °21 

q( 0 q( 0 
°21 '-’ll 


Noting from Fig. 4.1 that Ji = — T 2 , one obtains 


c(f) 

11 1 - r| e i2Af h ’ 

(f) _ (1 - r|) e iX < e < 
21 ! _ p2 e i2V f ’ 


r 2 (e i2Af£f - 1) 


(4.11) 

(4.12) 


where 


Z c + Zf 


(4.13) 


Finally, the composite ^-parameter matrix for a chafed cable may be written using the S'-parameter 
matrices derived above for nominal and faulty cable segments, with a switch to and from the corresponding 
T-parameter matrices. The result is 


S = 


c(f) „i2A c d f 

On e 

5 (f) e iA«(^ f ) 


c(0 P i\c(e-e,) 
°21 e 

^(f) e i 2 \ c (e-e,-d,) 


(4.14) 


4.2.3 System S-parameters 

The S'-parameter framework detailed in [78] is now applied to obtain the overall system model for the vector 
network analyzer (vna) based measurement setup shown in Fig. 4.2. Our ultimate goal here is to validate 
the fault model of §4.2.2; therefore, the two inner conductors of the STP cable are assumed shorted together 

t Alternatively, Tf may be computed as a single integral around a contour C that encloses both conductors. 
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Figure 4.2: VNA hardware model. 


at the ends in an effort to preferentially excite the common mode. The input signal is the source voltage vs , 
presented with a port impedance Z$. The load impedance Z\_ represents an arbitrary, possibly frequency- 
dependent, load connection, including the possibility of connecting the other end of the cable back into the 
VNA, in which case Z\_ = Zs- The output signal is the voltage indicated in the figure, which differs from 
the actual voltage at that point by an internal measurement delay t m and a calibration gain G. 

For modeling purposes, the VNA hardware measurement setup shown in Fig. 4.2 may be represented by 
the abstract S-parameter block shown in Fig. 4.3, with S given by (4.14). The general S-parameters for the 
composite system, including impedance mismatches, measurement delay, and gain, are given by [78] 


S = 


Q e iwt M 

i + s n r s — s 22 ri — (S11S22 — s 12 s 21 )r 5 r L 
Sn + r$ — r s r L s 22 — (SuS 22 — Si 2 s 2 i)ri 

s 21 (i + r s )(i + r L ) s 22 - r L 


5i 2 (i-r s )(i-r L ) 

IsTlSu + {SnS 22 — s 12 s 21 )r 5 


(4.15) 


where S mn are the elements of the common-mode S-matrix given in (4.14) for the possibly chafed cable. 

The elements of the matrix S may be measured directly by a VNA; equations (4.1) through (4.15) relate 
these measurements to the underlying cable and fault parameters that one may wish to infer from laboratory 
data. With only a look-up table for Tf computed ahead of time, these closed-form equations allow one to 
compute the system S'-parameters within ~ 20 /zs per measurement frequency oj m in MATLAB™ on a modern 
desktop PC. 


4.3 Bayesian parameter inference 

For the validation of the theoretical model against experimental measurements, one must adopt a method 
for estimating the effective parameters of the particular cable under test, as well as those of the measurement 


a 1 


a 1 


Zs 


Zc 
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Figure 4.3: S-parameter block diagram of the measurement setup with impedance mismatch. 
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environment. In this section, a specific probabilistic framework is presented for inferring the unknown or 
uncertain parameters 9 of the cable and the experimental setup from S'-parameter measurements at M 
discrete frequencies. In the most general case of all four system .S'-parameters simultaneously available, the 
noisy measurement model takes the form 

y = g(9,u) + v, (4.16) 

where u> £ R M is the vector of measurement frequencies, y £ C 4M is the vector of S'-parameter measurements, 
v £ C 4M is a complex vector of additive random noise, and g £ C 4M is a nonlinear function of 9 whose 
form has been elucidated in §4.2. The explicit w dependence on g will henceforth be dropped, but is implied 
throughout. 

Assuming statistically independent, identically distributed Gaussian noise for both the real and the 
imaginary parts of the measurements, the likelihood distribution for the measured data may then be written 
as I 

p(y\9, s) = (2ns 2 )~ 4M / 2s , (4.17) 

where the variance s 2 characterizes the uncertainty due to both modeling error and measurement noise. 

It is reasonable to assume that one would have access to manufacturer’s specifications, baseline or cal- 
ibration data, etc., from which one can gain some prior knowledge of the cable and experimental setup 
parameters. Let this knowledge be represented in the form of a joint prior distribution p(9) over the param- 
eters to be estimated. Additionally, the combined standard deviation s of modeling error and measurement 
noise is treated as initially unknown, and assumed to follow an exponential distribution p(s) = (l/s)e _s / s , 
for s > 0, with mean s = 0.1. The net prior distribution is then p(9 , s) = p{9) (1/s) e~ s ^ s . 

Ultimately, our interest lies in estimating the parameter values given the measured data. This is achieved 
first by specifying the posterior probability density function through Bayes’ theorem, 


P{d,s\y) 


P{y\9,s)p(9,s) 
f p(y\9', s') p(9' , s') A9' ds' ’ 


(4.18) 


and then finding the maximum-posterior estimate (9 * , s*) by solving the optimization problem 


maximize p(9, s|y). 

Furthermore, the shape of the marginal posterior distribution 

p{0\y) = J P{0,s\y) ds 


(4.19) 


(4.20) 


around 9* indicates how confident we are in this estimate. 

Unfortunately, finding an analytical solution to problem (4.19), or for that matter the marginal integral 
(4.20), is intractable. However, Markov-chain Monte Carlo (mcmc) methods enable one to draw N samples 
from the posterior distribution, 

9 n ,s n ~p(9, s\y), n = 1,2,..., AT, (4-21) 

with only the ability to evaluate the numerator in (4.18). This set of N samples then enables the estimation 
of the desired quantities. In particular, 

( 9*,s *) ~ argmax 0ii>Sii p(y\9 n , s n ) p(9 n , s n ). (4.22) 

Furthermore, by themselves, the 9 n samples are distributed according to the marginal posterior p(9[y), and 
thus their spread may be taken as a measure of our uncertainty about this estimate. More generally, the law 
of large numbers guarantees that 


£/(*») = [ f(0)p(9\y)d9 = E[f (9)]. (4.23) 

Thus, the samples can be used to estimate the expected value of almost any function of 9. Standard examples 
are the mean f{9) = 6 and the variance f(9) = (9 — E [9]) 2 . In addition, for skewed distributions with separate 

tSince y is a 4M-dimensional complex-valued vector, this is the joint distribution of 8 M real-valued random variables. 
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mean and mode, it is useful to represent the single- variate marginal uncertainty as a 100(1 — a)% posterior 
interval. This is defined as the interval [a,b] such that 100(a/2)% of the samples fall below a and above b 
[79, p. 38], 

There are many different MCMC-based algorithms one can implement to achieve the above sampling. The 
results presented in §4.4 were obtained using a relatively new method called nested sampling [80, 81]. 


4.4 Experimental results 

The modeling and Bayesian inference process outlined above is now demonstrated by a few examples using 
a 1-m long section of impedance-controlled aircraft-grade STP cable, whose geometric specifications are 
provided in Table 2.1, with SMA connectors attached to each end, and with the inner conductors shorted 
together at the source and the load ends in order to enable excitation of primarily the common mode. 
The S'-parameter measurements were made using a Rhode & Schwarz ZVB8 VNA (300 kHz - 8 GHz, 4 
ports). For all examples, the VNA was calibrated using a lower-grade kit (50 H, ZV-Z132, male), and unless 
noted otherwise, the S'-parameters were measured from 2 MHz to 8 GHz, in increments of 2 MHz, with an 
iF-bandwidtli (or measurement bandwidth) setting of 1 kHz and 0 dBm input power. 

The results presented below demonstrate that the twist-corrected analytical model can fit experimental 
data within 2% maximum error out to 1 GHz. Beyond this frequency, the (unmodeled) effects due to 
the connectors as well as cable imperfections start to dominate the data. This is a useful practical result, 
since most field-deployable instrumentation would operate within this frequency range for a variety of other 
reasons. We note in passing that the frequency-domain measurements reported here may also be used to 
synthesize the time-domain reflectometry (tdr) and time-domain transmissometry (tdt) responses of the 
cable, as we will show subsequently. 

4.4.1 Nominal cable examples 

Due to a variety of reasons including manufacturing variability in cable construction and calibration drifts of 
measurement equipment, it is advisable to infer the actual parameters of a given cable and experimental setup 
directly from laboratory measurements. Table 4.1 lists the prior information available for these experiments; 
the notation A f + denotes the nonnegative Gaussian distribution. The joint prior distribution p{9) of the 
vector of parameters to be inferred is given by the product of the marginal prior distributions indicated in 
the table. 


Table 4.1: Parameter prior information. 


Parameter 

Distribution 

Description 

£d 

A7+ (2.30, 2.0) 

effective relative dielectric permittivity 

<Jw 

A7+(6,0.2) x 10 7 S/m 

wire conductivity 

0’s 

A7+(0.5,2) x 10 7 S/m 

shield conductivity 

Od 

Af + (2, 2) x 10~ 15 S/m 

effective dielectric conductivity 

t 

A7+(l, 0.005) m 

cable length 


A/” + (50, 2) n 

source impedance 


J\f+ (0.1, 0.05) ns 

measurement time delay 

G 

S(G — 1) 

system gain (assumed fixed) 


Note that, in place of e and a used in §4.2.1, two effective parameters £d and Od are being introduced 
here to account for the inhomogeneous nature of the dielectric medium between the wires and the shield of a 
real impedance-controlled cabled If conductors 1 and 2 in Fig. 2.1 are assumed to be coated with insulation 

§ Since air and foam are nonmagnetic, the dielectric permeability remains //.(, . 
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of permittivity e = 2.3 eo and radius q (such that they are just touching each other and the shield), and the 
remaining space inside the shield assumed filled with air, a rough estimate for the effective relative dielectric 
permittivity would be 

e d ^ 1.57877. (4.24) 


Example 1 

We start by considering a simple field-deployable measurement setup where the input of the cable is attached 
to the VNA, and the other end of the cable is left open-ended in air. This is a more practical interrogation setup 
because it does not require attaching both ends of the cable to the VNA, or any specialized termination at the 
other end of the cable. (Generally, there are no impedance standards for aircraft connectors.) Furthermore, 
this measurement setup is the frequency-domain equivalent of the traditional tdr setup commonly used in 
practice. In this example, the VNA measures the system ,SVi parameter, which includes the signal reflected 
from the open far end of the cable, and thus provides more information regarding the unknown cable 
parameters than would be the case for a matched termination. Accordingly, the VNA measurements are 
modeled by 5n with 71 = 1 in (4.15). 

Table 4.2 shows the results of the model parameter inference for this example. The estimates and 
uncertainty information were computed from the approximate mode and 95% posterior interval of 2,000 
samples from the posterior distribution, all obtained using the nested sampling algorithm mentioned in §4.3. 
A cross-section of the samples, corresponding to the marginal posterior distribution p(e d, er w |j/), is plotted in 
Fig. 4.4. Perhaps most notable is the inference of the effective relative dielectric permittivity £d- The prior 
for this parameter corresponds to a homogeneous insulation of relative permittivity 2.3 (see Table 4.1), but 
the probabilistic inference algorithm was able to glean from data a more realistic value closer to the crude 
estimate given in (4.24). 

In Fig. 4.5, the measured data are plotted against the best (ie., most probable) fitting model for the 
system Su parameter. While only the VNA measurements out to 1 GHz were used to obtain the fit, the 
comparison between the measurements and model extrapolated out to 2 GHz is also shown. We also explored 
fitting all the data out to 2 GHz, and that resulted in increased error for the frequencies below 1 GHz. Under 
the controlled laboratory test conditions, it appears the residual error between the best fitting model and 
the measured data is more deterministic than random. This is caused by a variety of factors not considered 
in the model, such as frequency-dependent VNA miscalibration and geometry variation along the length of 
the cable. Furthermore, this error was treated as random noise in the probability model as an effective way 
to handle its presence. While the random noise component may be small in the lab, we do not expect this 
to hold in any practical environment. The modeling and parameter inference approach we have taken is 
expected to be robust under both types of uncertainty likely to be encountered in the field. 


Table 4.2: Example 1. Parameter estimates and 95% posterior intervals. 


Parameter 

Mode 

95% Interval 

Units 

fd 

1.5824 

[1.560,1.603] 


CT W 

6.0232 

[5.648,6.369] 

10 7 S/m 

0’s 

0.4974 

[0.479,0.520] 

10 7 S/m 

Cd 

2.000 

[0.156,5.699] 

10- 15 S/m 


50.8202 

[50.48,51.18] 

n 

t 

0.9997 

[0.993,1.007] 

m 

tu 

0.0776 

[0.077,0.078] 

ns 

s 

0.0063 

[0.006,0.007] 



52 




Example 1: MCMC estimation of £d and 
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Figure 4.4: Samples from the marginal posterior distribution p(Ed, a w \y), which show the uncertainty asso- 
ciated with the simultaneous estimation of these parameters given the measured data. 


Example 1: Modeled and measured Su 



Figure 4.5: Comparison between the measured and modeled 5?[(S'n)] (the plot for the imaginary part looks 
similar). Only the data out to 1 GHz were used in the inference process to obtain the fit, below this frequency, 
the maximum absolute error is less than 0.0213. The data beyond 1 GHz are provided to show the onset of 
unmodeled effects above this frequency. 
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Example 2 


The more ideal case is now considered where measurements of the system ,S' | i and S 21 parameters are 
available simultaneously. This will tell us how much improvement we can expect over the non-icleal measure- 
ment setup discussed in Example 1. In addition, the additional S 21 data will provide practical information 
regarding signal transmission through the cable, which is important for fault detection. 

The same table and set of figures shown for example 1 are produced for this example (also based on 2, 000 
samples from the posterior distribution). Table 4.3 and Figure 4.6 show that the estimated parameters and 
the associated accuracy were largely unchanged. This makes sense from an energy/information conservation 
perspective. In the first example, all the energy was reflected by the open end and measured at the input 
port. In this example, most of the signal energy is transmitted into the output port 2, and is measured as 
S 21 , while the rest is reflected back to the input port 1 and is measured as 5n. This example further shows 
that more measurements — there are twice as many measurement points in example 2 — do not necessarily 
translate to greater accuracy, a fact captured by our Bayesian inference approach. 

The most striking difference between the two examples occurs with respect to how well the model fits 
the measured data, shown in Figs. 4.7 and 4.8. The maximum absolute error over all measurements for 
both Su and S ' 21 was 0.018, and when the model is extrapolated beyond 1 GHz, the residual error grows 
much less rapidly with frequency. This might be expected because the cable is now properly terminated 
into the calibrated port 2 of the VNA, which was not the case in example 1. Another interesting difference 
is that the estimated port impedance parameter Z% differs from the estimate obtained for the first example 
by a statistically significant margin. This could be caused by a variety of reasons such as the non-ideal 
termination used for the first example or a possible drift in the calibrated impedance between the two sets 
of tests. 

Finally, Fig. 4.9 summarizes the practical results for the modeling work by showing the attenuation of 
a propagating mode as a function of frequency corresponding to the best model parameter fit given the 
measured data. The solid blue curve shows the attenuation constant a(co) = 3[A c (w)] for the model best 
matching the measured data out to 1 GHz; this is the attenuation one can expect for the given cable under 
perfectly matched conditions. The dashed red curve in the figure is the theoretical — In 1 5*2 1 1 (which is 
directly comparable to a), calculated by fitting the model to S 21 data measured out to 8 GHz in the lab, 
the latter plotted as the solid green curve. Since the S'-parameter model and data include a significant 
impedance mismatch between the source and the cable, higher attenuation is observed in an oscillatory 
fashion for frequencies greater than about 10 MHz, as expected in accordance with Fig. 3.6. The broad 
spectral feature peaking around 4 GHz in the lab data is believed to be a connector artifact; not surprisingly, 
the model does not fit this feature. 


Table 4.3: Example 2. Parameter estimates ± one standard deviation 


Parameter 

Mode 

95 % Interval 

Units 

£d 

1.5443 

[1.521,1.566] 


(Tw 

6.0201 

[5.634,6.374] 

10 7 S/m 

0's 

0.5003 

[0.480,0.526] 

10 7 S/m 

o-d 

2.000 

[0.169,5.742] 

10 -15 S/m 

Zs 

51.5423 

[51.19,51.93] 

n 

i 

0.9995 

[0.993,1.007] 

m 


0.0785 

[0.078,0.079] 

ns 

s 

0.0047 

[0.004,0.005] 
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Example 2: MCMC estimation of £d and er w 
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Figure 4.6: Samples from the marginal posterior distribution p ( ea , cr w |y), which show the uncertainty asso- 
ciated with the simultaneous estimation of these parameters given the measured data. 


Example 2: Modeled and measured Sn 



Figure 4.7: Comparison between the measured and modeled SR[(5n)] (the plot for the imaginary part looks 
similar). Only the data out to 1 GHz were used in the inference process to obtain the fit. 
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Example 2: Modeled and measured S 21 



Figure 4.8: Comparison between the measured and modeled real part of S 21 (the plot for the imaginary part 
looks similar). Only the data out to 1 GHz were used in the inference process to obtain the fit. 



Figure 4.9: Comparison of modeled and measured attenuation constants for STP cable (see the text). The 
cyan curve is the attenuation constant of theoretical RG-58 coaxial cable, provided for comparison. 
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4.4.2 Faulty cable example 

In this example, a 3 mm x 6 mm fault centered at 60° was added to the cable at 0.33 m ( i.e ., Wf = 3 mm, 
If = 6 mm, <f> f = 7t/3, and df = 0.33 m), which was then interrogated with the same measurement setup as 
in §4.4.1. The nominal cable parameters are assumed to have the values estimated in the previous section, 
and the goal is to validate the chafe model against laboratory data. 

In Fig. 4.10, the difference between the nominal and the chafed cable is plotted in the frequency domain, 
where we again focus on the real parts of Sn and SVi ■ The plot clearly shows increasing fault detectability 
with frequency, especially with respect to the Sn parameter, which is roughly equivalent to the TDR response 
of the chafed cable. Since the twist-corrected cable model matches the lab results out to about 1 GHz, and 
the chafe produces little measurable effect below 200 MHz, it is concluded that 200 MHz - 1 GHz is the most 
useful range of frequencies for model-based inference of unknown cable or fault parameters from experimental 
data. 




/ (GHz) 


Figure 4.10: Frequency-domain comparison between nominal and chafed cable models. 

In order to aid intuition, the time-domain transmitted and reflected responses to the input signal shown 
in Fig. 4.11 are synthesized and plotted in Fig. 4.12. 
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Time Domain Frequency Domain 



Figure 4.11: Input interrogation signal used for the time-domain analysis. 


tdr response 



tdt response 



Figure 4.12: Time-domain comparison between nominal and chafed models with lab measurements from a 
chafed STP cable. The plots show the net measured voltage response to the input pulse shown above in 
Fig. 4.11. 
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Specifically, the applied input source voltage was 


v s (t) = exp 


/ t — 1 ns \ 2 
^ 0.2 ns ) ’ 


Vs(f) = fft[u s (t)], 


(4.25) 

(4.26) 


where t = m At, for m = 0, 1, . . . , M — 1, with M = 8001 and At = 1/(2 M x 10 6 ). The net reflected and 
transmitted signals are then synthesized, from both modeled and measured data, according to the equations 


v R (t) = if ft 
vj{t) = ifft 


(1 + S u)V s 
<§21 Vs 


(4.27) 

(4.28) 


These eciuations follow directly from the definition of S mn , as described in [78]. The time-domain plots show 
the contribution from the chafe to be roughly the derivative of the input signal (as expected) and quite small 
at ~ 5 mV for the tdr response, and ~ 1 mV for the tdt response. 
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Chapter 5 


Conclusion 


5.1 Summary of main results 

From a modeling perspective, perhaps the most important point to stress is that the filament model is more 
than capable of reproducing the eigenvalues and the mode fields of both SUP and STP cables with sufficiently 
high accuracy for the purpose of fault detection. In fact, none of the other approaches that were attempted 
proved amenable to analytical treatment to nearly the same extent as the filament model. Some of the more 
esoteric features that are characteristic of helical and other periodic structures — significant difference in 
the frequency dependence of common- and differential-mode eigenvalues, presence of “extraordinary” bands 
(referred to as stop-bands in open periodic structures) with unusual dispersion and loss properties, etc . — are 
not observed in these cables over the spectral range of interest. It should be emphasized that this is due 
entirely to the very slow rate of twisting found in typical aerospace-grade STP cables, and is not a limitation 
or defect of the filament model. 

The essential consequence of this slow rate of twisting is a small correction to the eigenvalues and mode 
fields of SUP cable, which can be analyzed exactly and in great detail. The correction is important for 
propagation over long segments of STP cable; however, in analyzing local interactions — connection to source 
and load, scattering from a hole in the shield, etc . — the simpler SUP cable model may be used instead. 

The simulations have shown that the eigenvalues simulated from the simplified theoretical model are close 
enough to the eigenvalues of the true impedance controlled cable to be realistic for practical application. 
The fields generated by the historical tape model are not distributed anywhere near as uniformly as those 
in the wire model, and probably should not be used for reflectometry applications. 

The fields within a fault, which when scaled by the effective area of the fault are proportional to the 
expected reflections, are heavily dependent not only on distance from source but also upon the angular 
location. For a fault at a fixed distance from source, the fields were shown to vary by a factor of at least two 
due to angular variation alone. The calculated electrostatic 2-d capacitances varied by a few picofarads with 
respect to angle. The capacitances varied by almost 30 pF with respect to hole size ranging from 0.5 mm to 
2.5 mm. The ability to detect a fault can be thought of as being proportional to this change in capacitance. 
Most reflectometry systems will require at least a few picofarads of change to reliably detect a fault. 

Classical trade-offs still exist. For example, reflectometry that employs a fast rise-time edge pulse will 
better be able to discern the edges of a small hole fault the higher the frequency becomes, with the down- 
side that the loss with respect to cable length increases with frequency. The loss characteristics for shielded 
twisted pair are not really that different from coaxial cable such as RG-58. The difference is in the additional 
ambiguity that the angular position of the hole fault adds to the quantification of detectability. 

Modeling of a chafed segment of cable as a characteristic-impedance jump in an otherwise uniform line 
is a proven method that leads to a fast look-up table approach for fault detection [82]. Likewise, Bayesian 
inference of unknown cable and fault parameters using nested sampling was successfully implemented by the 
present authors before [78]. This approach met with great success here as well. Specifically, optimal fits to 
experimental data have demonstrated that the SUP cable model indeed provides an acceptable fit to data 
from STP cables up to 1 GHz. Not surprisingly, the slight analytical correction for twisting does not have 
a significant effect on the fits, but it does seem to reduce the bias on underlying parameter estimates. It is 
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therefore judged that other (unmodeled) elements in the system ( e.g ., cable imperfections, connectors, VNA 
calibration drifts, etc.) are mainly responsible for the residual. 

It is observed that the simple impedance-discontinuity model of a chafe is able to replicate the tell-tale 
“differentiation” effect of a small hole on a propagating pulse. It is also established empirically that average- 
size chafes do not produce a sensible reflection or transmission signal below 200 MHz, while unmodeled 
effects due to cable and connector imperfections start to become prominent above 1 GHz. It is therefore 
suggested that the theoretical cable and fault models derived here are suitable for chafe fault detection within 
a spectral window of roughly 0.2 - 1 GHz. 


5.2 Areas of future research 

Various open questions remain: 

• first-principles modeling, verification, and validation of chafe faults, particularly when both cornmon- 
and differential-mode signals are present on the cable; 

• analytical modeling and/or numerical simulation of distributed (random) scattering from cable imper- 
fections {e.g., rough shield surface, nonuniform twist rate, etc.)', 

• development of closed-form or look-up-table S'-paranreter models for typical connectors; 

• physics-based trade-space study of fault detectability in practical settings; 

• (more) realistic noise models for more robust fault detection algorithm design; 

• implementation of fast deterministic optimization algorithms for on-line fault detection; 
among others. Some of these items are currently in the investigation pipeline. 
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